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I. Introduction 

This annual report contains investigations on the properties of 
nonideal plasmas, which were carried through in the period from 1 November 
1978 to 1 November 1979, under ONk Contract N00014-79-0073. Progress 
was made in the evaluation of the electrical conductivity and the 
thermodynamic functions of nonideal plasmas. Research which was pub¬ 
lished in this period is not included. 

Chap ter II gives a dimensional analysis of the possible electrical 
conductivities of nonideal classical and quantum plasmas, with and 

without thermal effects. Quantum effects are considered since the 

21 -3 

electrons are noticeably degenerate for densities n > 10 cm and 
temperatures T * 10^°K. In strongly nonideal plasmas, y - 7.e ? n ^/KT >> 1, 
the thermal energy is negligible compared with the energy of the 
Coulomb microfields, i.e. such plasmas behave like zero-temperature 
systems. 

Chapter III presents a kinetic theory of the electrical conductivity 

of nonideal plasmas, when the electrons have i) a Maxwell and ii) a 

Fermi distribution of velocities. This work is based on the quantum 

mechanical scattering cross section for an effective, shielded 

Coulomb potential, which is applicable to intermediate nonideal conditions, 

0.1 < y t 10. For moderate degeneracy, n < n, the conductivity 

depends on a Coulomb logarithm, but not for complete degeneracy, 

3/ ] 5 i/ 

n>>h, where n = 2(2nmKT/h 2 ) ■- = A. 828 x!0 T ■ . 

Chapter IV is concerned with a statistical theory of Lite free 
energy of nonideal gaseous (0 y < 1) and quasi-liquid (1 < y *" y ) 
plasmas, where y c is the critical interaction parameter of solid 
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plasma. The physical model used considers short and long range Coulomb 
interactions through quasi-lattice interactions and collective electron 
and ion waves. The degeneracy effects on the free energy are discussed 
for large interaction parametersy. 

Chapter V derives the Hamilton function and the canonical field 
variables for a many-component plasma continuum with longitudinal 
Coulomb field interactions. By Fourier analyzing the canonical 
fields for random fluctuations the Hamilton function can be used to 
calculate the free interaction energy of nonideal plasmas. This 
method has,however, not yet lead to concrete statistical applications 
because of the mathematical difficulties associated with the 
evaluation of the complex integrals in the multi-dimensional phase 
space of the Fourier amplitudes of the canonical fields. 

Chapt er V I gives a simple application of the theoretical approach 
in Chapter V concerned with the evaluation of the distribution function 
of the velocity fluctuations in a neutral one-component gas. The 
theoretical distribution function is shown to be in agreement with 
the observations in turbulent gases. 

Ch apter V I1 is an appendix, in which an unrelated subject is 
discussed, namely the propagation of stress relaxation waves and their 
use for signal propagation and system detection in water. 

The main results obtained are summarized in the abstracts of 
Chapters II - VII. Other research concerned with calculations of 
the electrical conductivity 1) using quantum-field theoretical methods 
of solid state physics and 2) considering the effects of the fluctua¬ 
tions of the Coulomb microfields on the current transport could not be 
written up in time, and will be communicated separetelv. The same holds 
for tin; Ph.D. thesis of Mr. A.H. Khalfnoui on nonideal plasmas. 
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II. Dimensional analysis of electrical conductivity 

OF NON IDEAL CLASSICAL AND QUANTUM PLASMAS 


Abstract 


By moans of dimensional analysis, novel formulas for the electrical 
conductivity of nonideal i) classical and ii) quantum plasmas are derived 
based on the axioms of Dupre"! In the general case of a nonideal plasma 
with partially degenerate electrons, the conductivity is of the form a * 
C q > T ^ y B en ^ /ra ^ 2 , where C q is a dimensionless constant, A and B 
are powers, and y T * e 2 n ^/KT and Yq * e 2 n ^/ (fi 2 /ran ^) are the reduced 
(7!=1) nonideality parameters of the classical and quantum plasma, 
respectively (e, m, n are the charge, mass, and density of the electrons, 
KT is the thermal energy, and it is Planck's constant). The known 


conductivities are obtained as special cases of this conductivity 

3/ 1/ 

r ormnin, e.g., the conductivity a ■ C (KT) 2 /m /2 e 2 of the ideal, classical 

o 

plasma, or t lie conductivity a * C^(me 6 /fi 3 )n ^/KT of the high-temperature 


met a1. 
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I 


INTRODUCTION 


In recent years, numerous measurements of the electrical conductivity 

1 — 8 ) 

of non-ideal plasmas have been reported,- which were produced by oven 

heating ^(medium pressures) and shock wave heating —^ (high pressures) of 
alkali vapors and noble gases, with typical pressures ranging from 10* 
to 10~* bars. In spite of the availability of an approximate kinetic 
equation for nonideal plasmas which considers spatial and temporal 
correlations in the collision operator, satisfactory theoretical explana¬ 
tions of the experimental conductivity data on nonideal plasmas are 
missing to date. The degree of the nonideality of a fully ionized plasma 

is measured in terms of the (dimensionless) interaction parameter y, which 

, 1 /, 

represents the ratio of average Coulomb interaction (Ze z n J ) and thermal 
(KT) energies (n = electron density, Z = ion charge number, e = elementary 
charge), 

y = Ze 2 n ' 'KT . 


-3 -1 

In cgs-unit, y = 1.670 x 10 Zn J T . The conductivity theory of ideal, 

12 ) 

fully ionized plasmas— agrees with the experimental data only if y << 1. 
For moderately, 0.1 < y < 1, nonideal plasmas, the ideal conductivity 
theory yields much too large conductivity values-—- (KT) ^ 2 /m ^ 2 e 2 Z fn A, 
where l\ = {l+(D/p Q ) 2 =D/p q for D >> p Q , p^ = Ze ? /2KT is the impact 
parameter for 90° deflections, and D is the maximum impact parameter. 

The electric shielding length of Debye is related to the interaction 
parameter y by 

D = [Z/4n(l+Z)] ^ 2 y ^ 2 n 


The ideal conductivity theory breaks dc'm at higher electron densities 
because tiie Debye radius D loses its physical meaning as an electric 
shielding length and upper impact parameter when the number of electrons 









in the Debye sphere, = 4nD 3 n/3, is no longer large compared with one. 
and y are related by 

SI D = (4tt/3) [Z/4tt (1+Z) ] ^ y ^. 

For typical nonideal conditions, n > 10 2 °cm 3 and T = 10 4 °K, the Debye 
_8 

radius is D <10 cm, i.e., is smaller than the atomic diameter, which 
shows that the electric shielding concept is not applicable to proper 
nonideal plasmas, y > 1. Another reason for the inapplicability of the 
ideal conductivity theory to proper nonideal plasmas is its assumption 
of successive, small binary interactions, whereas in reality a conduction 
electron experiences many-body interactions for y 1. 

In the following, we apply dimensional theory to the derivation of 
new formulas for the electrical conductivity of (non-relativistic) ideal 
and nonideal, classical and quantum plasmas. In the most general case 
of an electron plasma, the electrical conductivity a is a power function 
of the characteristic plasma parameters, a ~ e P m CJ ti r (KT) S n t (tf is Planck's 
constant). As special cases, the conductivity formulas for the ideal, 
classical, fully ionized plasma and the partially degenerate, solid metal 
are obtained. The derived formulas for nonideal classical and quantum 
plasmas indicate the dependence of o on y, which can be compared with 
the experimental observations. 
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THEORETICAL FOUNDATIONS 

In a system of reference in which magnetic fields are absent, a linear 

4 -4 

electric current response j =oE exists, provLded that the generating 
electric field E is sufficiently weak. For any gaseous, liquid, or solid 
medium, the electrical conductivity o = |j|/|F,j is given by 

a = (nc 2 /m)x (1) 

where e is the charge, m is the mass, n is the density, and t is the 
(average) momentum relaxation time of the current carriers. Because of 
the large ion mass m^ >> m, the main current carriers m a plasma are 
the electrons. Eq. (1) holds for any perturbed Maxwell or Fermi distribution 
of the electrons. 

• 13) „ 

Dimensional analysis is based on the axioms of Dupre. - By axiom 
1), absolute numerical equality of quantities a and b may exist only when 
the quantities are similar qualitatively. That is, a general relation 
may exist between two quantities a and b only when the two quantities 
have the same dimension. By axiom 2), the ratio of the magnitudes of 
two like quantities a and b is independent of the units used in their 
measurement, provided that the same units are used for evaluating each. 

In general, any measurable quantity o(the secondary quantity) 
can be expressed in terms of tiiose appropriate quantities a^, i = 1,2,... M 
(the primary quantities), which affect the magnitude of o. The general 
relationship between the magnitude of the secondary quantity a and the 
magnitudes of the primary quantities a^ is a function of the M arguments 
of the form 

o = f(a j, a 2 , a 3 ,...a M ) 


( 2 ) 
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Application of the axioms 1) and 2) to Eq. (2) demonstrates that the 


functional relation fCa^ is the power function- 
Nj N2 N3 N . 

0 = C a *a z a„ ... a, M 
0 1 2 3 :l 


1A) 


(3) 


C is a dimensionless coefficient, which depends on the nature of the 
0 

physical quantity 0 , and can only be determined by means of a detailed 
physical model. In many cases, the order-of-magnitude of is one, 

C '1. 

o 

In the most general nonrelativistic case of a thermal quantum 
plasma, of which the classical thermal plasma is a special case, the 
secondary conductivity quantity 0 depends on the dimensional primary 

quantities a[ = e(clectron charge), a ; > = m(electron mass), a 3 = ffU’lanck 1 s 
constant), a [( = n(electron density), and a,. = KT(t!iermal energy). 

The dimensionless constant C Q is in general a function of the dimensionles 
parameters p/p^ of the plasma, 

C = C fZ., m/m.,...p/p.,...) . (A) 

0 01 1 1 

E.g.,Z is the ratio of the magnitudes of the ion and electron charges, 
m/nu is the electron to ion mass ratio, etc. The electrical conductivity 
0 and its primary quantities have the following dimensions P(L = dimension 
of length, T = dimension of time, M = dimension of mass) and units U in 
the cgs-system: 


Via] 

= T -1 

♦ 

U[a] 

-1 

= sec 

V[e] 

= 

* 

U[e] 

3 /? 

= cm ' gr 

V(m ] 

= M 

9 

U[m] 

= gr 

W 

= ml 2 t _1 

9 

m] 

= gr cm 2 

P[n] 

= IT 3 

9 

U[n] 

= cm 

P[K") 

9 -2 

= ml 2 t 

9 

U [ KT ] 

- gr cm 2 


( 5 ) 
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It is recognized that the secondary and primary quantities depend only 
on three basic dimensions which are independent, namely L, T, and M. !or 
this reason, dimensional analysis provides at most three independent 
equations for the determination of the powers , i = 1,2,3,...M in Eq. (3). 
That is, at most three powers can be calculated while at least M-3 powers 
have to be determined by comparison with experiments or by physical 
arguments. 

Two of the primary quantities are variables, namely n and KT 

J ry Q 

(K =• 1.380 x 10 gr cm"sec / K), whereas the remaining three primary 

-io Vo Vo -l 

quantities are elementary constants (e = 4.803 x 10 cm ^ gr sec , 
m = 9.109 x 10 8 gr, tl = 1.054 x 10 2? gr ctn 2 sec ). 

The fundamental equation (3) is applied below to the determination 
of the electrical conductivity of (nonrelativistic) non ideal, classical 
and quantum plasmas, in which the electrons are responsible for the 
electric current transport. To illustrate the results, they will be 
expressed in terms of the interaction parameter y for Z = 1 and the order- 


of-magnitude of fundamental electron energies: 


y t = E c /E ? = ne 2 /KT , 

E j = KT 
E c = e 2 n V3 

_2 / 

= fi 2 /mn 3 

is the order-of-magnitude of the quantum potential energy Q = 

, , .. 1 L 

-(Tv /2m)V (, z /() of the electron in the plasma. E and E^, are the 


(b) 

(7) 

( 8 ) 
(9) 


Coulomb interaction and thermal energies. 


CLASSICAL T=0 PLASMA 


In a classical. (fi>-o) electron plasma, in which the thermal energy KT 
is negligible compared with the Coulomb interaction energy, the conductivity 
depends on the dimensional parameters e, m, and n. By Eq . (3), 


Nj N 2 

a = C e *m n 

o 


Hence, 


T' 1 - , 


( 10 ) 


( 11 ) 


(3/2)N, - IN, = 0, (1/2)N! + N 2 = 0, -Ni = -l. 

These are three independent equations, which determine the powers 
Nj, N 0 , and N 3 uniquely, 


( 12 ) 


N, = 1, N = -V. , 


N i = * ■ 


(13) 


( 14 ) 

}h 


By Eqs. (10) and (i 3) , the conductivity of the zero-temperature, 

classical plasma is 

o = C en ^ /m ^ 
a 

Eq. (14) indicates that a = C a <i^/(4Tr) ^ increases proportional to n /2 , 
1/ 

wiiere w = (.4Trne : ’/m) ■’ is the plasma frequency. This result was first 
P 

derived by Euneman— and latpr by Hamberger and Friedman —' for an 

electrostatically turbulent T = 0 plasma by means of semi-quantitative 

1 / 1 / 

physical arguments, which give C^ = ;r ' (in^/m) for - I . 


i 
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CLASSICAL T>0 PLASMA 

In a classical (h>o) , thermal (KT>0) electron plasma, the con- 
ductuvity depends on the dimensional parameters o, in, n, and KT. 

By Fq. (3) 


o = C e Nl m N;y n N ^ (KT) N| ' 

0 


Hence t 


T~ 1 - (lVV‘)% N ^*3 ( „|/ T - : ' ) S| ’ , 


l .e. , 


JN] - 3N 3 + 2N 4 = 0, ^-Nj + N 2 + Ni, = 0, - N, - 2N, ( = -1 . 

These are three independent equations, which determine three of the 
powers in terms of the fourth. 


Nj = 1 - 2N, N 2 = - }• , N 3 = 


_ 1 _ 1 


N, N,, = N 


(15) 

(16) 

07) 

four 

(18) 


According to Eqs. (15) and (18), the conductivity of the classical, 
thermal plasma is 

« - C o o'- 2 ".,-* „ 5 - 5" (KT)" . (19) 

Collecting of powers of N reveals the dimensionless group contained in 
Rq. (19), and condenses the conductivity formula to 

a = C ( (KT/e 2 n ^ 3 )^ en^/m^ 2 . (20) 

For an ideal (classical) plasma, o = (nc 2 /m)t cannot depend on n 
since x ^ n 1 for binary e-i collisions. Hence, a and N are for an 
ideal T • 0 plasma 

= r ti rrt h h 


a = C (KT) 7 /e 2 m ■ , N = 3/2 , 


( 21 ) 


] 2 ) 

in agreement with kinetic theory—- , which shows that C 
1 / — 1 

3/4 (2xr) ^ZCnA ~ 10 . For N = 0, Rq. (20) reduces to o of the classical 


T = 0 plasma, Rq. (14). 







For a nonideal (classical) plasma with v-body interactions, we 
have t “ n 1 V and a - (ne^/m)r * n ? V , with v > 2, i.e. N = 3v - 9/2 -• 1/2. 


Thus, we find for the conductivity of nonideal, classical plasmas 

„ -N V? , V ? ... 

o = en / /m ' , N ■ 3/2 , 

where Y, r = E../E t ' le non ideality parameter defined in Eqs. (6) - 
l 1 


( 22 ) 


(S). F.q. (22) expresses the important result that the conductivity of 

-N ty> 

a nonideal, classical plasma, o = /2rt 2 , decreases proportional to 

-N 

Y^- witii increasing nonideality since N > 3/2. The exact value of N 
can bi determined by comparison with experimental data. 





1, 


T = 0 QUANTUM plasma 

In a completely degenerate electron plasma, F^. •< Eq, the conductivity 
depends on the dimensional parameters e, m, n, and fi, but not on KT. By 
Eq. (3), 


„ N, No N,.N 4 
o = C e 1 m ^ n • ft H 


Hence, 




l . e. , 

jN[ - 3N 3 + 2N 4 = 0, jNj + N 2 + N 4 = 0, -N, - N 4 = -1. 

These are three independent equations, which determine three of the 
four powers in terms of the fourth, 

N, = 1-2N, N 2 = -- - N, N 3 = i + In, N 4 - 2N. 

Substitution of Eq. (26) into (23) yields as conductivity of the 
completely degenerate electron plasma 


a = C 




v* 


For N = 3/2, Eq. (27) leads to the conductivity of a solid metal 


at T = 0, 


= C tf^n/e ? m^ , 


N = 3/2 , 


17) -I 

where • C <* Z . For N=0, Eq. (27) reduces to a of the classical 


T=0 plasma, Eq. (14). 

For complete degeneracy, E^, is negligible compared with E^ ~ -mVp , 
1 / 

where Up - (tf/m)n * is the Fermi velocity. For this reason, Eq. (27) 
is rewritten in the form 

o = C (j en ^/m ^ (29 


where 


Yq = K ( ./Kq = e ? n ^/ (1v /mn ^) 


(30) 
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is the nonideality parameter of the completely degenerate electron plasma 
for Z = 1. Since N = 3/2 for the T=0 metal, it is to be expected that 
N>0 for the completely degenerate electron plasma, i.e. its conductivity 
decreases with increasing nonideality y , Eq. (29). This formula is 
useful for the interpretation of conductivity data of completely 
degenerate electron plasmas, with N as adjustable parameter. 
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T -0 QUANTUM PLASMA 


In a partially degenerate electron plasma, E,j, £ E^, the conductivity 

depends on the dimensional parameters e, m, n, KT, and fi. By Eq. (3), 

„ N, N ' N 3 N 4 i .N 5 

o^C^e '-n J (Kl) 4 fi J . (31) 

Hence, 

T~ * = (L^M 1 / 2 T - 1 ) Nl M N 2 L~ ,N 3 (ML 2 T _? ) N ''(ML ? T _) ) N,i , (32) 

i . e . , 


-N! 

- 3N 3 + 

2N 4 

+ 

2 Ni; = 0 

9 



+ N 2 + 

n 4 

+ 

N 5 = 0 

9 

(33) 

-N, 

- 2N„ - 

n 5 

= 

-1 

. 



These are three independent equations, which permit to express three 
of the five powers in terms of the remaining powers, 

Nj = 1 - 2A - 2B , N 2 = - i - B, N 3 = 5 - |a + |b, 

N 4 A, N 5 = 2B . (34) 

Combining of Eq. (31) with (34) gives for the conductivity of the 
T>0 quantum plasma 


■=• w (^) B 


(33) 


whore A and B are powers which can not be determined by dimensional 
reasoning. 

For A = -1 and B = -3/2, Eq. (35) yields the conductivity of solid 


metals at temperatures T>0, 
me 6 n 1/s 

a = C —-— , A = -1, B = -3/2, 


0 *3 


(36) 


KT 


where— a Z ^. Eq. (36) expresses the 1/T - lav^^of the 
metallic conductivity at "high temperatures". 

Eq. (35) contains the dimensionless groups and Yq» which permit 
to rewrite the conductivity formula as 

A -B V 2/ V 2 

en /m . 


0 = C o Y T Y Q 


(37) 



Since A = -1 for T>0 metals [Eq. (36)] and B = 3/2 for T=0 metals[Eq. (28)] 
one can speculate that A > -1 and B>0 for nonideal, partially degenerate* 
plasmas. A and B can readily be determined by means of conductivity 
measurements for nonideal quantum plasmas. The theoretical determination 
of the nonlinear dependence of o on y^ and Yq in Eq. (37) from a 
physical model of more than two-body interactions is left to future 
research. 

In conclusion, it is noted that we have derived new conductivity 
formulas for nonideal classical T>0 plasmas [Eq. (20) or (22)], completely 
degenerate plasmas [Eq. (27) or (29)], and nonideal T>0 quantum plasmas 
[Eq. (35) or (37)]. These formulas can be used to interprete conductivity 
measurements on nonideal plasmas. Once the still undetermined powers 
N, A, and B are known empirically, it should be possible to develop a 
conductivity theory for nonideal plasmas which provides an explanation 
of the y and y. dependence from first principles. 
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III, Electrical conductivity of nonideal plasma 


ABSTRACT 

The electrical conductivity of fully ionized, moderately nonideal. plasmas 

2 173 

with interaction parameters 0.1 < y S 1, where y ■ Ze n /KT is the ratio 

of Coulomb and thermal energies, is calculated for displaced Maxwell and 

Fermi electron distributions, respectively. The electrons are scattered by 

-1 -1/3 

an effective Coulomb potential 4>( r ) * Zer exp(-r/<$), with 6 * (3n/4nZ) ' 

the mean ion distance, which considers binary (0< r $ 6) and many-body 
(<$ < r < ») interactions. It is shown that the resulting conductivity formula 
is applicable to densities up to four orders of magnitude higher than those 
of the ideal conductivity theory, which breaks down at higher densities 
because the Debye radius loses its physical meaning as a shielding length and 
upper impact parameter. In the limit of complete degeneracy, the conductivity 
formula reduces to that of a solid metal. 
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.1-3) 


lations^ give an electrical conductivity of the form a ~ (KT) ^ / ^/m^ / ^e‘'Z('ttA 


INTRODUCTION 

The theory of the electrical conductivity of fully ionized plasmas' J/ based 
on the Boltzmann equation, the Fokker-Planek equation (derived by expanding the 
binary collision integral for Lhe small, successive velocity changes of 
Coulomb scattering), or the Lenard-Ba1escu equation (taking into account the 
dielectric properties of the medium) is in agreement with the experimental data 
for rarefied high-temperaLure plasmas, y<<1. The interaction parameter is 
defined as the ratio of (average) Coulomb interaction (Ze^n^^) and thermal (id) 

energies (n is the electron density and Z the ion charge number), 

2 1/3 -3 1/3 

y = Zcn ' /KT = 1.670 x 10 Zn /T 

in cgs-units which will be used throughout. The conventional transport calcu- 

i form a 

2 1/2 

for classical ideal plasmas, where A = [1 + (D/p n ) ] ~ l)/p for D>^p o . D is 

2 1/2 

the maximum impact parameter (Debye length), D = fKT/4ire (1 + Z)n] , and p 

2 

Is the average impact parameter for 90° deflections (Landau length), p = Ze /2KT. 

The condition, A>>1 or £nA~1.0^ is satisfied only for not too low temperatures 

4) 

T and not too high densities n. Conductivity formulas with this Coulomb 
logarithm break/down for large interaction parameters y and densities n, since 
Lhe Debye radius 

l>- |Z/4.(1 W.)l 1 / V' / V l/3 

“8 

becomes smaller than the atomic dimension 10 cm and, thus, completely loses 
its physical meaning as an electric shielding length and maximum impact parameter. 
E.g., tor T = 10^ °K, y>-10^and D<10 ^cm if n>10^cm Moderately nonideal 
plasmas with y~l are readily generated through shock wave compression and exhibit 
conductivities of the order u 10^ - 10^mho/cm^ which are much smaller than 
those which would be obtained by applying the conductivity formula for ideal 
plasmas in the nonideal regime. 
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Although there are some bulk measurements of the eleetric.il conductivity 

5 — 8 } 

of nonideal cesium and noble gas plasmas available , theoretical explanations 

of these results are still missing. The momentum and ncrgy transport in 

9) 

weakly nonideal plasmas,y « 1, was treated by Wilhelm by means of an 
experimentally shielded Coulomb potential, which permits Lo consider not only 
short-range binary (r < D) but also long-range many-hody (r D) interactions. This 
interaction model ;as used shortly afterwards by Rogov for the calculation 
of the conductivity of weakly, nonideal argon and xenon plasmas with Debye 
shielding. 

For moderately nonideal plasmas, 0.1 <y S 1, various phenomenological approaches 

have been used to extend the conductivity formula of ideal plasmas, e.g., Goldbach 

et al multiply the Debye length D with a free parameter x(p)which is chosen 

to match the experimental data, i.e. to compensate for the too rapid decrease of 

D with pressure. A kinetic equation has been proposed for nonideal plasmas 

12 ) 

by Klimontovich , which considers spatial correlations and temporal retardation 
in the collision integrals. This equation has not yet lead to transport 
coefficients because of the mathematical difficulties associated with its solution. 

In the following, the momentum relaxation time and the electrical conductivity 
of (i) classical and (ii) degenerate plasmas is calculated for intermediate non¬ 
ideal conditions, 0.1 < > S 1. For this region of interaction, the concept of 


Debye shielding already breaks down since the number of particles in the Debye 

3 

sphere 4irD /3 is no longer large compared with one for y > 0.1. This difficulty 
can not be remidied by replacing D with the quantum mechanical shielding length 


(even when thermal effects are included)’ 


D p = (na o /k F ) 


a -Sly, k, 
o mo 2 1' 8,1 


which is of the same order as D in most high pressure plasmas, e.g. D ~ 10 cm 

r 

for n = lO^ cm From the definition of the mean particle distance, it is clear 


that the mean ion distance 6 ~ n separates the region in which an electron 
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experiences few - body encounters (r < () from the region in which an electron 

experiences many - body interactions (r > 5) in a nonidea] plasmas, as long as 
—8 24 — 3 

5 > 10 cm (n^ < 10“ cm ). Thus, the mean ion distance evolves naturally 
as the characteristic interaction distance for nonideal plasmas, for which 
Debye and Fermi shielding fail. 

We calculate first the electrical conductivity of plasmas with (i) Maxwell 

and (ii) Fermi distributions of the electrons, when all ions have the same 

charge number. Then, we generalize the conductivity formulas for plasmas 

with several ion components of charge number Z^, i = 1, 2, ... . The electrons 

are assumed to be scattered by the exponentially decaying Coulomb potential 

p = Zer ^ exp(-r/<5) which takes many - body interactions at distances r > 6 

into account. The considerations are applicable only to moderately nonideal 

24 

conditions, 0.1 < y < 1 , up to densities n << 10 cm. 

24 -3 

For solid state densities and larger densities, n > 10 cm , we have 
-8 

6 < 10 cm and the chosen interaction potential t( r ) is no longer valid. For 

24 -3 

densities n >10 cm , the Coulomb field of the ions is changed by polari¬ 
zation of the bound electron cloud, so that the free electrons are scattered 
as in solids by the ions. For this reason, plasmas in the solid phase are 
not treated herein. 
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PHYSICAL FOUNDATIONS 

The electrical conductivity o of ,-my gaseous, liquid, or solid medium, in 

which the electrical current transport is due to electrons, is proportional 

to the electron density n and the relaxation time t of the average momentum 

<mv of the electrons (m is the electron mass and e > 0 is the elementary 
o J 

charge) 

2 

o = (ne /nt) x ( 1 ) 

hq.(l) reduces the calculation of a to the evaluation of i. The relaxation time 
r is determined by Lhe scattering potential and the (classical or quantum 
statistical) kinetics of the electron gas in the electric field. 

In proper nonideal plasmas, the region 0 < r < 6 of binary and few - body 
collisons and the region 6 < r < « of many - body interactions are bounded by 
the mean ion radius, 

(4n5 3 /3)n. =1, 6 = (3Z/4tm) l/3 , (2) 

since electric Debye shielding exists only for weakly nonideal conditions, y << 1. 
For this reason, the effective Coulomb potential of Z times charged ions is in 
plasmas of intermediate nonideality 

<1>(r) = Zer ^ exp(-r/ 6 ), 0 < r < <® , 0.1 < y < 1. (3) 

Eq.(3) is no longer applicable to plasmas with densities comparable to solids. It 

contains the binary and few - body collisions at distances 0 < r < 5 and the 

many - body interactions at distances 6 < r < 

The differential cross section o(0,g) for the scattering (g -*• g *) of 

14) 

electrons by the potential (3) is in the center of mass system 

o(0,g) = (Ze 2 /2m) 2 /[g 2 sin 2 (0/2) + u 2 ] 2 (4) 

where 

0 = Mg,g*), g = v c - Vj_, 8* = v* - v £ * (5) 

and 


u = Ti/2m4 


( 6 ) 







Tile electron and ion velocities oefore and alter llie interaction are designated 
-y -> 

by v . and v* . , respectively. The speed u corresponds to a de Broglie 

G ) 1 G t l 

wave length of the order \ - 6. For u >• 0 or ■ - Kg.(4) reduces to the 

Rutherford cross section. 

The scattering cross section u(0,g) is strictly valid only in the Born 
approximation^^. Contrary to what one might expect in general lor the latter, 
Eq.(4) describes in good approximation the scattering in the exponentialiy 
decaying potential (3) because of the peculiarity of the Coulomb interaction. 

The Coulomb interalion (J> - 1/r has the unique property that the Born approximation 
and tlie exact wave mechanical approach give the same scattering cross section^ 
(identical with the Rutherford formula). In the region 0 < r < 6, the 
interaction potential (1) is practically Couiombic, and thus the Born approxi¬ 
mation gives the correct solution. In the region 6 < r < », the interaction 
potential (1) is effectively screened, i.e., the Born approximation gives the 
correct solution because ^(r) is small. 1n the transition zone r ‘ n, the 
Born approximation bolds fairly well for reasons of continuity. 

The relaxation time i is obtained by evaluation of the collision integrals 

-► 

for tlie electron momentum mv for the (i) classical and (i i) degenerate plasma, 

respectively. Both in the cases of classical and Fermi statistics, the particle 
-y >■ 

velocities v . and v* . before and after the interaction are interrelated by the 
tj* i- g , i 

conservation equations for momentum and energy. 




CONDUCTIVITY OF CLASSICAL i’LASMA 


According to kinetic theory, tin- average momentum density n m(-'v ■ 

e 1 

exchanged per unit Lime between electrons and ions, interacting, with the Coulomb 
potential tl), is given by the collision integral lor mv , which determines 

U 

the momentum relaxation time t , 

-> • > 

“Illll(<v > - <V.>)/l “ 

o 1 ' 

m /.../ V |f (V*) 1 . (v*) - f (V ) f. ( v.) ] ,; .(',g) d:,; dv dv. . (7) 

J 1 e e e 1 l e e i i e i 

The scattering cross section c (>, g) is given in L- • j. C 4) .and the sol iti angle element 

is d'-? = sin 0 di) do. In response to an applied electric fit Id E, Lite electrons 
t ► >• 

and ions dr-'ft with velocities <v and <v. > so that their distribution 1 suctions 

e i 


are displaced Maxwellians, 

f ,('■ ,) = n (m ,/2-tKT ) ^~ oxp[-~ m (v - <v . •)‘7KT ) 

S S S S S iL S !■> , t s 


■ , i. 


(«) 


E<(. (8) represents a 5-moment-approximnt ion to the nonecpii 1 ibrium solution of Lite 
Boltzmann equation.' 1 *' The perturbations of f (v ) due to viscous stresses and 

S ii 

heat flows are neglected in Eq.(8), since they yield only corrections of 
higher order to the conductivity. 

9) 

The collision integral (7) is integrated by standard methods for subsonic 
drift velocities, j <v > - <v > | < (2KT/ra)* with the usual approximat ions 


(m . = m m. / (m + m .) 
e L e l e l 


m : in, T = m I (T /m ) + (T./m.)] T 
e e S t s' e e J . .. 


T; . 


For supersonic drift velocities, a linear response j = oL between current d. nsity 

9) 

j and electric field E does no longer exist. The resulting relaxation time is 


m 


given by: 

-1 8 .. 


"(2KT/iim) 1/2 u. (> 


^(Ze 2 / KT ) 2 L 


L, (A ') 


(h) 

( 10 ) 

OD 


wlieri 


I’KT/mu 2 - <Hiii/'ir ) ( m/ 1Z) V, | 


( 12 ) 




is t ho exponential integral of order one (T = 0.477 ... - Killer's constant). ' 


The latter satisfies the inequalities, for x 0, 


+ 2/x) < e X E L ( x ) - t'»U + 1/x), (1 + x) ' ,.'i j(x) < X ( 

According !y, E<|. (1 1) gives formally lor sm.il] and 1 urge :> reunion tu = A K 

L ■- tn A , A > 1 , ( 

r. - a , a « i ( 


Combining of Eqs. (9) - (12) with ]•; |. (1) yields tlx- desind electric 
conductivity of the classical plasma of iniormed i.ile nonideality, (J. I < c < 1 : 

o = 3(K.T)- J/2 /2(2tim) 1/2 e 2 El. ( 

where 

- •’ _ > / a 

h = f'n[8mti ~ (4ttn/3Z) KT |, A >• I, ( 

-2 -2/3 

1. = Binlf (4ir n /3Z) ' KT , \ « \. ( 

by Eqs. (15) - (16).It is instructive te rewrite A in terms of the thermal and 
quantum potential energies of an electron, 

A - 8E T /E Q , E t = KT, E = ir/mh 2 . ( 


The conductivity formula (17) 
plasma^ mainly through the term 
logarithm, L = Cn A , for A ■ • 1, i 
which the plasma is nondegenerate, 

A = 3.482 x K) 11 (n/Z) _J/3 T 
The corresponding argument 2KT 


differs from the conductivity of the ideal 
L. The latter has the form of a Coulomb 
.e. for all densities n and temperatures T 
E.j, E , Kq. (20) . Numerically, 

3 1-3) 

D/Z.' - of L lie ideal Ciuilomh logarithm , 


for 


( 








Table 1 compares A of the nonideal plasma and of the ideal plasma for large 

densities n and the typical temperature i ; 10 ,J K. It is seen that hi A 

i) of 

the ideal plasma is unacceptably small for densities n > I (J ' 8 cm ^, whereas 

, 22-1 
t II ' ol the nonideal plasma has reasonable values up to densities n ■: 10 cm , 

if T = 10* °K. The conductivity formula (17) holds, therefore, for densities n 

up to A orders of magnitude higher than the conduct ivitv formula of the ideal 

plasma. l’.q.(17) is not applicable to n - T r. j-.utvi for which A [, L.e. 

E *< F. , Eq. (20) , since in this case the i I • t i mi;; would he degenerate. 


nlcm J ] 


TAB LE 1 : A and A versus n f. r 1' - 10 °K and 7. = 1. 


3.482 x 10 3 1.616 x 10 2 0.750 x 10’ 0.348 x 10° 

1.035 x L0 1 1.035 x 10° 1.035 x 10 _1 1.035 x 10~ 2 


The conductivity formula (17) becomes in cgs - units or practical units 


(y x 10^ see 3 = 1 mho citt'b, 


• = l.'394 x 10 8 T 3 ^ 2 /Z hi A [sec ‘ ] = 1.549 x 10 4 T 3,/2 /Z hi A [mho cm ] ] (23) 
where A Is given in Eq.(21). Accordingly, if 1 = !0 + °K and Z = 1, o = 1.899 x 10 
mho enf 1 for n = 10 J8 cm” 3 and o - 2.097 x 10 1 mho cm -1 for n = 10 2 ° cm 3 . 









CONDUCTIVITY 01' QUANTUM IM.ASMA 


The electrons in a plasm.-, become degenei ‘i .v if their thermal DeBrogl. ie 


vave lenj’.th is larger than the mean electron I i si .nice, i.e. al densities (h=2iiti) 


n • 2 (Turn KT 


/h 2 )3/2 = 4ji28 x 10 15 T 3/ "\ 


2l ..„-3 


ll.g., for T = 1O* 4 K, degeneracy requires n ; x \(r tin . In view of their 


arge mass m. 


m, the ions can be treated as classical. The momentum relaxation 


time i of the degenerate electron gas is determined by the quantum statistical 

-► 19 ) 

collision integral for niv , 

e 

-> • -V 

- nm(<v -• - <v.a)/i = 


v {f (v*)f.(v*)[l - f (v )] - f (v ) f. (v.) [ I - Vf 'v*)l 
J J e e err 2 e e e e t l 2 e e 


x g o(0,g) dW dv^ dv^ ( 

where the scattering cross section o^.g) between electrons and ions is given by 
Eq. (4). The solutions to the velocity distributions are the displaced 
Maxwellian (8) for the ions (s = i) and tiie 5-moment Fermi approximation for the 
electrons, 

f (v ) = 2(m/h) ! tl + exp[Jan(v - v •)" - u|/KT i 1 ( 

e e 2 e e 

The chemical potential p = p(n, T) is determined by the integral functional 


n = J f „(c ; p) dc . 
e e e 


Again, a linear response j = oE exists for small drift velocities <v 

e, i 

or weak electric fields E. integration of Eg.(25) yields, after standard 

approximations, for the relaxation time of the degenerate electron gas: 

-1 _ 8.2 1/2 e 4 Z 2 nj L Q __ 

3 m (RT) 3/2 R(n,T) 


where 




. _ 16 m 4 an.-2/3 n 

\) ~ 3 ^2 ( ~3Z ) Kf a(n ’° 




2 7 


and 


Q(n,T) = |(1 + 2 2 " + ...), n <• n (1') , (30) 

2- n 



R( n,T) = 

ii 

(1 + 2 2 ~ + ...). n n(T), 

n 


(31) 

but 





Q(n,T) 

3 3n 1/2 n. 
* ' n 

m \i + - 5 r(- 3 4- ?)■*'’•* - .1. 

^ n 

it ' • n(T) , 

(32) 

R(n,T) 

1 3n 1/2 n 

0 * /. .: ' 

, q 2 . 1/2 -4/3 

[ 1 + ~ (--;-- -:) 

n - • n(T) . 

(33) 


Equations (30) - (31) and Eqs. (32) - (33) result from expansions of the Eermi 
distribution (26) in the collision integral (25) for: densities n « n(T) and 
n ->'■ n(T) , where 


2 3/2 

n (T) = 2 (2-nm KT/h") J/ 


(34) 


is tire critical density which separates the degenerate and non-degenerate regimes. 
These series are based on expansions of the normalization integral of Eq.(26), 
which gives the chemical potential m explicitly as a function of n and T, 

-£= = t»l{" [1 + 2~ 3/2 (") 1 + 4 " 3" J/2 )(") 2 +...], n ' n(T), (35) 

Kin n ^ n 


P_ _ , 3 ii 
KT ' A 


0 ... 2 , 1/2 . A 1/2 H/ „ 

n. 2/3 . _n , 3n _ n.--A/3_* .Jh__ n -8/3 

1 12 1 4 „ 80 A ' „ 


n - n(T) (36) 


n ^ n 

Combining of the conductivity formula in Eq.(l) with the relaxation time of 

Eq(27) yields for the electrical conductivity of the degenerate electron plasma 
of intermediate nonideality, 0.1 < Y :> 1: 


3(KT) 3/ 2 R (n,T) 
8(2m) 1/2 e 2 Z L 

where I. is given by Eq.(29), 


(37) 


In the limiting cases of large and small values of 


I., = 0 (n ,T) | , 

» > i 

h () = - 2 A0(n,T), 


A 0 
0 


(38) 


(39) 
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since A = (2/3)AQ by comparison of Eqs.(29) and (12). 

For n/n -> 0, Eqs.(37) and (38) reduce to the. classical conductivity, Eqs. 
(1/) and (18), since R(n,T) > 4//n and 0(n,T) <• 3/2 for n/n 0 by Eqs. (30) 

- (31). On tlie other hand, Eqs. (37) and (3')) p.ive in the limit of coniplete 
degeneracy, n/n ■>■ 


0 "9 ~2 2 

2 inn e Z I, 


where 


, =3 3 L/3Z_ 

Q 40' , 2ir 3 


l q ;<1 


by Eqs. (32) and (33). Euqations (40) - (41) combine to the conductivity lormula 


.15tt 2 .2 a 1/3 h^n 

.,6 C 3 ’ 2 2 _5/3 

2 m e Z 


Equation (42) agrees with the expression for the conductivity of a low temperature 

, 20 ) 

metal.— 

For numerical evaluations, Eq.(42) is stated in cgs units (sec ) and 


practical units (mho cm ), 

_ /. /.no ia’^ n -5/3 


LB „ —3/3 


j - 4.498 x 10 1' n [sec ) = 4.998 x 10 Z a [mho cm |. 

Accordingly, a 4 x 10^ mho cm -1 for n = 10^ cm' J and u 5 x |o mho cm ' 
for n = lO*" 4 cm 3 ? if z = 1. 
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GENERALIZATION 

Nonideal plasmas exhibit frequently not onl\ a high degree of 

single ionization but also multiple ionization, due to lowering of 

the ionization energies by the internal Coulomb fields, and over 

lapping of the atomic wave functions at sufficiently high pressures. 

In an electrically neutral plasma with N species of ions (i) of 

charge Z.e and density n^, the electron density n and entire ion 

density n(i) are related by 

N N 

n = l ^| n i > n(i) = l n. . (44) 

i=l i=l 

The characteristic interaction radius 6 of the Coulomb field of each 
ion is within the many-component plasma 

6 = f 4im(i)/ . (45) 

Since the probabilities for interaction between the electrons 


and ions of type i = i, 2,...N are additive, the momentum relaxation 
time of the electrons is 


where t. is the relaxation time of the momentum exchange of the 

l ' 

electrons with ions of type i. 


Classical Plasma. By Eqs. (1), (46), and (9)-(10), the conductivity 


of the many-component plasma with Maxwell electron distribution is 
in the Lorentz approximation 


2 (2irm) • e 


3 (KT) 
-1/— N- 

I '« V 


(nj/n)Zj •!. 
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I, is given by Eq. (11) in general, and the approximations (15) and 
(16) in the limits A >> 1 and A << 1, respectively, where now 

A = 8KT/(Ti 2 /m6 2 ) . (48) 

Quantum Plasma . By Eqs. (1), (46), and (27), the conductivity 
of the many-component plasma with Fermi electron distribution is in 


the Lorentz approximation 


3/ 

3(KT) 2 R(n,T) 


!; N 

8(2m) '• e 2 T (n./n)Z.*L 

i=i 1 1 Q 


Lq is given by Eq. (18) in general, and the approximations (38) 


and (39) in the limits A^ >-■ 1 and A^ 


.1 , respectively, where 


A q = (16/3)KTO(n,T)/(-ri 2 /m6’) . (50) 

For a plasma with only one ion component, we have N = 1, 

7,^ = Z, n./n = Z 1 , and Eqs. (47) and (49) reduce to the previous 
formulas (17) and (37), respectively. In case of the many-compnent 


plasma, not only the electron density n but also all ion densities n^ 
have to be known for the evaluation of the conductivity. 
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IV. Statistical thermodynamics of nonideal plasma 


Abstract 

A quantum statistical theory of the free energy of a nonideal electron-ion 
plasma is developed for arbitrary interaction parameters 0 < y < y 

c 

1 /. 

(y = Ze'n J /KT is the ratioof mean Coulomb interaction and thermal 

energies), which takes into account the energy eigenvalues of (i) 

the thermal translational particle motions, (ii) the random collective 

electron and ion motions, and (iii) the static Coulomb interaction energy 

of the electrons and ions in their oscillatory equilibrium positions. 

From this physical model, the interaction part of the free energy is 

derived, which consists of a quasi-lattice energy depending on the 

interaction parameter y, and the free energies of the quanitized electron 

and ion oscillations (long range interactions). Depending on the degree 

of ordering, the Madelung "constant" of the plasma is c t(a) = a for y>>l, 

Vo 

u(y) = ct for y > 1, and a(y) <r y " for y << 1, where a - 1 is a constant. 

The free energy of the high-frequency plasmonr. (electron oscillations) is 
shown to be very small for y > 1, whereas the free energy of the low- 
frequency plasmons (ion oscillations) is shown to be significant for 
y > 1, i.e. for proper nonideal conditions. For weakly nonideal plasmas, 
y " 1, both the electron arid ion oscillations contribute to the free energy. 
Unis, novel results are obtained not only for proper nonideal (y > 1) 
but also for weakly nonideal (y << 1) plasmas. From the general formula 
lor f lie free interaction energy hr of the plasma for 0 < r ' ■ Y , simple 
analytical expressions are derived for AF in the limiting cases, y >> 1, 


1, and y ■-< 1. 





INTRODUCTION 


In the classical work of L)ebye and Htieckcl on electrolytes, the total 

Coulomb interaction energy is calculated from the continuum theoretical 

picture of every ion interacting with its surrounding space charge cloud. 

Using more sophisticated methods, similar results were obtained for 

1 2 

weakly nonideal plasmas (y<<l) by Mayer (cluster expansion), Ichikawa 

3 A 

(collective variable approach ), Vedenov and l.arkin (graphical density 
expansion), and Jackson^ (hydrodynamic continuum interaction model). 

Based on different methods and approximations, investigations of moderately 

(y > 1) and strongly (y>>l) nonideal plasmas were given by Berlin 

6 7 8 

and Montroll , The inter and Centry , Kcker and Kroell , Ebeling, Hoffman 

9 in 

an ' Kelbg , and Varobev, Norman and Vilinov , respectively. 

in spite of differences in the theoretical approaches, the. leading 

terms of the analytical results for proper nonideal plasmas (y > 1) 

give essentially the same 1 formula for the free plasma energy, AF/NKT = -ay + 

2 1/3 

blny + c, due to Coulomb interaction, where Y = Zc n /KT is the 

ratio of electron - ion interaction energy and thermal energy, and 

a, b, c are constants depending on the respective approximations 

and assumptions. The thermodynamic functions of strongly nonideal 

plasmas (Y'-->l)were also determined with the help of Monte Carlo and 

computer methods by Brush, Sahlin, and Teller^, Hansen^, I)eWitt^\ 

14 

and Theimer , respectively. Although computer methods provide limited 
physical insight, they are useful for checking the quantitative validity 
of analvtical theories . 


At sufficiently high electron densities, for which Y > 1, classical 
statist ical theories fail due to thermodynamic Instability^, which is 







inhibited by quantum mechanics. The classical plasma pressure would collapse for 
Y> 1 due to the negative electron-ion interaction energy, whereas in reality 
the pressure remains positive in a plasma due to the Fermi pressure (exclusion 
principle) of the electrons. For these reasons, we present herein a quantum- 
statistical theory for nonideal plasmas based on concepts similar to those 
used by Debye for solids^^. The application of this model to proper 
nonideal plasmas (y > 1) is justified since a plasma exhibits a quasi¬ 
crystalline structure for y > 0.1 before it undergoes a diffuse transition 

into a solid, metallic state at a critical value y . The roll of the 

c 

longitudinal phonos of the Debye theory is assumed by the quanta of the 
plasma oscillations (plasmons) in the case of the quasi-crystalline plasma. 

The theory is also applicable to weakly nonideal conditions, since the quasi¬ 
lattice energy reduces for weak ordering, y <■ 1, to the free interaction 

energy of weaklv nonideal plasmas. 

The theory to be presented takes into consideration (i) the energy 
eigenvalues of the random, collective electron and ion oscillations and 
(ii) the static Coulomb interaction energy (quasi-lattice energy) of the 
electrons and ions in their oscillatory equilibrium positions. Thus, 
all significant long and short range Coulomb interactions are considered. 

The results are applicable to arbitrary nonideal plasmas, 0 < y < Y c , where y c 
is the critical ordering parameter at which a phase transition into a 
solid metallic state occurs. 
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FtiYSlCAh : o;:N|)ATIi>t;S 


Subject, of the theoretical considerations are quasi-homogeneous high- 
pressure plasmas consisting, of electrons of charge -e and density n = N/V and 

ions of charge l/.e and density u/Z - K/ZV, with typical densities in the range 

2()-3 IL ( _j 3 4 0 

IlT cm < n •_ lir cm and temperatures of the order-of-magnitude T 10 - 10 K. 

2 

For these conditions, Lhe Debye radius D = [4iinc (1+Z)/KT] 11 is D = 6.901 x 

1 —B 

[T/n(1+Z) 2 •; 10 cm, i.e. , is smaller than the atomic dimension and the number 
of particles in the Debye sphere would be N 4tniD /3 < 1 for Ik 10 cm and 

>4 _3 

nv 10 cm . It is seen that the concept of Debye shielding completely breaks 

down, and statistical theories containing the Debye length as a characteristic 

parameter would be physically meaningless for high density plasmas. 

The nonideal behavior of plasmas is determined by the interaction parameter 

2 1/3 

'i , which is the ratio of the Coulomb interaction energy ~Ze n and thermal 


energy KT, 

y = Ze 2 n 1/3 /KT= 1.671 * 10“ 3 Zn' /3 /T. <D 

It follows that 0.5Z < y < 15Z for 10 2 ^ctn •_ n < 10 2+ cm 3 and T ~ 10^ °K. For 

Y = I, the nature of the plasma changes from a "thermally expanding" (y^l) to 
an "electrostatically contracting" (y>l) plasma. For y>l, the collapse of the 
plasma due to Coulomb attraction between electrons and ions is inhibited by Lhe 
Fermi pressure of the electrons, i.e. by the quantum mechanical exclusion prin¬ 
ciple. Thus, in the region 0 < y <■ y the plasma undergoes a diffuse transi- 

e 

tion from a nonideal classical plasma (y < 1) to a quasi-crystallinc. plasma 

(1 *■ y ■ y ), with an incomplete ordering comparable to that of a liquid, 
c 


An understanding of strongly nonideal plasmas has been attempted via the 

i • i (6-14) 

model ot discrete interacting particles in a dense gas - . lor the above 

reasons, however, it appears to be a more reasonable procedure to calculate 






tlif thermodynamic functions of proper non idea 1 |> J a: anas i i om Liu.- picture of 

col Uu't ive 1 ‘Ji‘i'lron ;r w l inn oscillations. In this app roach, the free inter¬ 
action energy is due to the static Coulomb interaction of the electrons and 
ions in their "equilibrium positions" (.‘ladeiung energy) and their oscillation 
energies about the equilibrium positions (plasmon energies). 

Since the plasma volume V contains ,\ electrons and N/k ions, there exist 
SN (high-frequency branch) and iN/Z (low I i equeiicy branch) characteristic fre¬ 
quent. ies 'X. ot longitudinal oscillations. Kacli plasma oscillator of frequency 
U. can only have the energy (n. + '. ) h. ,n -0,1,2,..., s> that the energy Fa i j 

l tit 

of a plasma state with n. plasmons ol frequency u. is 

l i 

K { i } “ ) n .diui. (2) 

I i} 1 1 

where {if designates the entire set ot given eigeninqaencies m. Ac¬ 
cordingly , the partition function Q of rlu- longitudinal plasma oscillations is 

t) = )’ e- K,i,/KT = II l/d-e^i 7 ^ 1 '). (3) 

{i > 1U 


Front Q, the thermodynamic functions such as the pressure, internal energy, 
entropy, etc., are derived in Lite usual way, e.g., the free energy of the 


plasmons is 


F - -KT In t) = K.T ^ ln(l-e , 


In the limit V*-" , t he discrete eigenf reipienc ies o>. are t'ep Laced by con¬ 
tinuous ones, ti>=id(k) , in accordance with the uispersion law for space charge 
waves of wave length = 2it/k, () y k < k. 

I. Fleetron Oscl llati ons . The high-frequency branch of the space charge 
waves i s cltio to longitudinal electron oscillations. Their i requency m is for 
classical (n '• n) and completely degenerate (n • n) electrons given by-*-^ 


ur [1 f(r /4 tt)ZY ' (kr ) ) 
]> e v 


1 




n • * n 


( 6 ) 


<7 


or - m 


11 + zsr (-^) 1/| <t> 7l zi'‘ l i k ' : u ) i. 


where 


. y, 

n = 2(2)imKT/lr ) 

1 / 

,,i = (Anne- /m) ‘ 

P 

r- =n-^ 

t? 


(7) 

(H) 

(9) 


are the critical electron density, the plasma frequency, and 

the mean electron distance (k = c /c of the electrons, and m is 

e p v 

their mass). Since k ~2ir/r (oscillations with \ < r are 

max e e 

physically inconceivable), the electron oscillations propagate, 
w = w(k) > w , in non ideal plasmas. 

2. Ion Osc ill a tions . The low-frequency branch of the space 
charge waves is essentially due to ion sound waves. Since the ions 

are presumed to be nondegenerate, the frequency of the ion oscillations 

, 17) 
is given by-- 


1 / 

u - A (k )(k .KT/K) ‘ k 


( 10 ) 


where 


' (k) 


1 + 


z(kq/< ;) 


VA 


1 + (k: /Air)ZY ' (kr )' 


( 11 ) 


<5 (k) = 1 , n • n , (1.2) 

is a correction factor of magnitue-of-order 1, which shows the 

influence of the electrons on the ion oscillations (M = mass, r. = 

c /c of the ions). 

P v 

In weakly nonideal plasmas, y • i, the electron sound waves are 


strongly damped for wave lengths X • D, due to trapping of the resonance 







electrons with thermal speeds comparable to the wave speed. For prope 
nonideal plasmas, y % 1, the number ol particles in the Dobvc splure 

, -H 

4nD J /3 is no longer large compared with one and D ■■ 10 cm is smaller 
than atomic size, so that thermal Landau damping is no longer feasable 
For this reason, electron oscillations should exist for wave Length 






VI 


ST AT I ST I f’.AI. TIIIT.MUDYNAM I ( S 

In the plasma iimli'r cons Ldernt ion, I lie e i oil run:; aiui ions interact 

through their longitudinal Coulomb fields (transverse electromagnetic 

2 

interactions are negligible for KT me ). 11 u • electrons (s - e) and 

ions (s = i) have tliermal velocities c and random collective mean 

mass velocities v due to their oscillatory wave motions about the 

s 

equilibrium positions, so that /;/ i f d"V = f), and • v -■ = 0 where 

I is the equilibrium velocity distribution of the- species s and 
s 

designates a spatial or temporal average. The resulting, Hamilton function 
with Coulomb interaction leads to a tree energy of the plasma of the 
form: 

1 ’ )' h 'J 0) + ii M + l ! s . (13) 

s=e,i s-o, i 


1 ^°^ is the ideal tree energy of tin non im eraet. i ng plasma components 
s. K is the Coulomb interaction energy ol the electrons and ion'; in their 
equilibrium positions. F^ . is the free energy of the electron and ion 
oscillations, i.e. of the high and low frequency plasmons, llq. (A). 

It. should he noted that bq. (’>) takes into, eons iderat ion all 
significant short-range and long-range Coulomb interactions by means ol 
the Made lung energy F.,, and the plasmon energies F . As is evident from 


17) 


the derivation of bqs. (5)-(6) and flO), in wiiieh terms of order m/M are 

neglected compared to 1, Fq. ( (1 ) coin ains the e-e, e-i, and i.-i Coulomb 

- 1/1 

interactions at distances A > n 

1. Free Energ y Fs °^. In high pressure plasmas, the electrons are 

153/2 -3 

part iallv degenerate for densitie: a where n = A.828 x 10 T |cm ], 


whereas the ions behave in general classically. 


for the free energy of the ideal electron gas 


8 ) 


Fermi statistics gives 




4<l 


i'^' 0 = - NKT U. 1/ .,(p/KT)/U 1/ (n/kT 


(14) 


whei\ 


U (i./KT) 


J_ C _x P dx _ 

(i’ + l) J ^x-m/KT f ! 


I/:’. •*/:? 


(13) 


and 


n = 2 (2irmKT/h“) 32 “ U^., (p/KT) 

19 ) 


(16) 


defines the Sommer f eld integrals , .and let <-rin i lies tie* eliemicnl 

potential p=p(n, T) of Lite eleetrons, respectively. The free 

energy of the translational degrees of freedom of I lie classical, 
18) 


ideal ion gas is 

F (o) = - (N/Z) KT ln[(2nMKT/h 2 ) 3/2 
\ 


(17) 


2. Quasi-Lattice Rnergy E^. The equilibrium positions of the 

electrons and ions, about which the electrostatic oscillations 

occur, form an electron "lattice" and an ion "lattice", with 

18 ) 

an incomplete ordering. By means of the Ewald method , one 
calculates the Coulomb ineraction energy of the qunsi-cubically 
centerid elect run-ion lat tices as 

K m = - re, NKT , <» = ct = 1.451 for Y -• l. (18) 

As the ordering of the plasma increases with y, a (y) is a weak 

function of y such that asymptotically r x = 1.541 for y 1. l r .q.(18) 

2 _ 

indicates that -E /N Ze /r Is of the order of the average 
M 

e-i interaction energy. For weak ordering, y ■ 1, it will he sho.Ti 

V, 

Y 


that <t 
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3. High-Frequency Contribution F . Since the number of Longitudinal 

modes with wave numbers between k and k + dk in volume V is V4 tt k 2 dk/(2 tt ) 2 , 

Eq. (4) gives for the free energy F of the high-frequency electron 

oscillations of energy Tim(k) 

k a 

F o /KT(V/2tt 2 ) = l!n{ .l-cxpl-ti'to(k)/KT) } k : 'dk (19) 

0 


where 


, , v, 

m(k) = u (1+a 2 k 2 ) '< , 


c z /ur = (k /47T)(Z/>)r z , n << n , 
m p e e 

(3/5) vJf/u 2 = (-~) ^(-v-) (Z/y)V , n >> n , 

1 p 2 Uti on e 


( 20 ) 

( 21 ) 

( 22 ) 


by Eqs. (5)-(6). The speed of sound and the Fermi speed Vp 


of the electrons are 

c = (k KT/m) 1/2 , 
m e 


v,, = 'ff(3Tr 9 n) ^/m. 

r 


The number of modes in (0, k fi ) and V equals the number 3N of 
degrees of freedom ot the electron gas, i.e.. 


(2r) ’V 


K 1 / 
r c 4irk 2 dk = 3N, k = (18irn) h 

J e 

0 


(23) 


(24) 


Integration of Eq. (19) by parts yields, under consideration of 
k ! KTV/hn 2 = 3NKT, for the free energy of the high-frequency plasmons: 


Tim 


3NKT bill - exp[ - 


JL 


KT 


( 1+a' k^ ) ^ | f - F ( 


Tim 


KT 


■ • lk 0 ) 


(25) 


who re 


lito hi.', ak i(/, , ;■, - ’/' , 

ak > ■ —' c <ak J " - 

KT KT (tfco /KT) (i+x 2 ) 1 

0 ° 


(261 


and 


li'.i /KT = ( 4 ;t) ^ (A /r ) (y/ i.) ^ 

p e e 


> = b/(mKT) ,/? , 

o 


( 27 ) 
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ak = k ^ (9it ^ ? /4) ^(Z/y) ^ , n <•- n , 
0 0 


ak £ = 2 ^ & 7i (3/2) (n/n) ^(Z/y) ^ , n n. 

By means of the successive substitutions, (f) x-sinln,, dx = cosh c, df 
and (ii) c = ("Rio^/KT) coshf,, dc » (Roj^/K' 1) slnh fd ( .,the integral (26) 
is transformed to 


F(e p’ ak e ) = (ak e' : p )-3 \ (f 2 _t? p) ( c ' - l) ‘ <l ’ 


Vo / _ 


where 

e = diw /KT, e = e cosh(arcsinh ak ). (3! 

p p ’ e p e 

Since the leading expression in Eq. (25) is the logarithmic term, it 

is sufficient to give for F(e^, a ^ e ) tl.e series approximation (Appendix), 


, 3 . -% , 


F(c ,ak )/2 ' 2 (ak ) 3 e 
P e e p 


l e mE p l („)(2e ) n m + n) y(^ + n, (c -e )m), e < 3e , (32) 

,= i n=n n P 2 e p e p 


where 


5 e -e 3 

, „ - + n e P 2 +n _ 

Y(- + n, (e -e )ra) = m / u e mu du 

2 e p J 


is the incomplete gamma function— . Since in general y/Z ~ 1 for 


Cp •- e^ < ^Ep, the expansion (32) is useful where simple approximate 
relations do not exist. 
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4. Low Frequency Contribution F With the number of modes in the 
interval dk at k and volume V given by V4rk 2 dk/( 7 it) 3 , Eq. (4) yields for the free 
energy of the low-frequency ion oscillations of energy Tfm(k) 


F./KT(V/2tt 2 ) 


In {1-exp [ --Rw (k ) /KT ] > k 2 dk 


where 

u>(k) = 6(k)c M k, 

1 / 

c M = (r.KT/M) /? . 

by Fqs. (10) and (12). The number of inodes in (0, k^) and V equals 

the number 3N/Z of degrees of freedom of the ion gas, i.e., 
k 

(2n)~ 3 V ( 4Tik ? dk = 3N/Z, k f = (18u ? n/Z) ^ 


Partial integration of Eq. (34) gives, under consideration of 
k? KTV/6ir 2 = 3(N/Z)KT, for the free energy of the low-frequency plasmons: 


F. = 3(N/Z)KT I fn{1-expf- 


6(k.)k.]) - G(k. 

i t j 


where 


(l(k ) = k - 3 f 1 U(k)4 k 6‘ (k)lk 3 dk 

1 KT 1 j e (fiCj i j(KT)6(k)k 


Since the dispersion factor <5(k) is a bounded function varying very 

1 L 

little with k such that 1< S(k) ^ (1+Z) ' for k r (0, k^), 6(k) can 

be approximated by an average value 6 , 

6 (k) = 6 1, n ^ n . 


Since in addition the logarithmic expression is the dominant term in Eq. (38), 
the integral (39) is approximated by 


G(e ) = e f 1 t'(e f -1)~ 1 dr 
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where 


e = "hc^Sk/KT, e L *-fic^k /KT. 


70) 


G(e^) has the serai-convergent series expansions, 

G( Ci ) = 3U— + 20 z i + - * -1 ’ e i << lf 
G( 6l ) = -jy + 0[e c i], e i >> 1. 

This completes the formal mathematical aspects of the theory, 
the physical implications of which require further elaboration. 


(42) 

(43) 

(44) 
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APPLICATIONS 

For applications of the theory to strongly, Intermediate, and weakly 

nonideal plasmas, it should be noted that the dimensionless parameters 

y/Z , 'flu) /KT, ak , and n/n occuring in Eq. (25) for the free energy F 
p e e 

of the high-frequency plasmons can not be varied independently. Since 
both y/Z and A__/r^ increase with increasing n and decreasing T, 

Tho /KT ~ (A g /r g ) (y/7.) ' 2 varies over a large n-T region similar to (y/Z) 1 , 

Eq. (7). Numerically, 

y/7 = 1.670x 10“ 3 n 1 / 3 /T, -fito /KT = 4.310x10 _7 n^ /T, n/n = 2.071xl0” 16 nT"^ 
ak > = 2.219 (y /Z) 12 , n << n 

ak g = 1.910(n/n) 73 (y/Z) % n>>ii. (45) 

E. g., for T=10 4 °K, y/Z£ 1 if n <’ 10 21 cm “ 3 and !•:,»/KT > 1 if 
n > 5xl0 2 °cm~ 3 . For T= 10 3 °K, y/Z < 1 if n i 10 K Vm" 3 , etc. Thus, 
for typical conditions of nonideal plasmas y/Z and Urn /KT are of the 
same order of magnitude. It is also recognized that in general 
n/n •■> 1 if y/Z >> 1, and n/n << 1 if y/Z<< 1. 


Tn Eq. (38) for the free energy F^ of the low frequency plasmons, 
only one characteristic parameter e^ occurs since 6(k)~ 4 ~ 1. By Eq. 
(42), this parameter is 



(18ir 2 ) 1/3 K i 1/2 5 -— = 2.158x10 5 Z ^ (•£}•) ^ ^ ^ K< 1 ( 46 > 

r. 

i 


where 


=-ft/(MKT)^ , v i = (n/Z)^ 3 . (47) 

Accordingly, for typical nonideal plasma conditions, it is << 1 
since > /r^* - ' 1 (classical ions) although in general A^/r^ > 
electrons) for y/7, > 1 or 1i«) /KT > 1. 



1 (degenerate 
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The deviation AF of the free energy of a nonideal plasma from 
ideality is by Eq. (13) due to the quasi-lattice energy F.^ and the plasmon 
energies ., 

AF = E w + Y F . (48) 

s=e, l 

Since the theory of electron oscillations^^has not yet been 
developed for arbitrary degrees of degeneracy (n > n), the contributions 
of the electron iscillations to AF in the cases n n and n ? n have to 
be estimated from the dispersion equations for n << n [F,q. (5)] and 
n >> n [F.q. (6)], respectively. Fortunately, it turns out that 
]F |<<|AF| for y/ z> l> so that quantitatively rolyable approximations 
for AF can be derived. 


1. Strongly Nonideal Plasmas. By Eq. (6) the spectrum w(k) of 

electron oscillations extends over a band Am - w above the plasma 

- < - P % -3 

frequency for y/Z » 1 since kr g - k^r ~ 1 and (n/n) ~ 1. Application 

of the mean value theorem for integrals to Eq. (25) shows that the free 


energy F^ of the high-frequency plasmons vanishes exnonentially for 
£ i. e. y / 7. -* oj : 


? /3NKT = (ln {i - exp[-c (l+a 2 k 2 ) ^ ] 
e y p e 


e (ak ) 
p e 

exp[Cp(l+x‘ 


—77- r x^l+x 2 ) ^dx^ v 0, c 

) ,/2 ]-l J / p 


0 - x - ak 

e 


Accordingly, | | /3NKT << 1 for £^>>1, i.e., y/7, •>!. On the other 

liand, the free energy of the low frequency plasmons is by Eq. (38) for 


nondegcuerate ions 
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F. = 3(N/Z)KT[£he, -(1/3)] = 

^ ^ V 

3(N/Z)KTU'hy + £n[ (18 tt 2 /Z 4 ) 1/3 ? ] - (1/3)}, e. <■<•' 1. (50) 

e ? /fi 1 

It is noted that y/Z >> 1 is compatible with e. = 1Tc^k^6/KT<< 1 as 
explained above. 

Equations (49) and (50) demonstrate that the contribution of the 
electron oscillations to the free energy is negligible in strongly 
nonideal plasmas, y/Z >> 1. In this limit, the nonideal part of the 
free energy is due to the quasi-lattice energy E^ and the ion 
oscillations, 

AF/NKT = - ay + (3/Z)6ly + (3/Z)£rt (f^A^) - 0/7.), y/Z » 1, (51) 

where 

v B = e 2 /h, 6 = (1 8 tt 2 Z~^) . (52) 

Note that tny depends on both n and T whereas fiiPc^/Vg depends only 

Q 1 I 

on T, where the Bohr speed is v fi = 2.118 x 10 cm/sec >> c M = (i^KT/M) ' 2 . 
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It is remarkable that the electron oscillations contribute little 

to the free energy compared to the ion oscillations for y/Z •> 1. This 

result holds even for moderately nonideal conditions, y/Z > 1. Thus, 

we disagree with the formula "F = ne + 3NKTf’lt (1'uj /KT>" stated without 

o o 

21 

derivation for nonideal plasmas by Norman and Starostin- , according 
to whom "all the vibrations have exactly the same frequency u 
near the plasma frequency < 0 ^". The derivation of this formula requires 
fiw(k)/KT << 1 for the electron oscillations, which implies y/Z << 1, 
but the latter inequality contradicts their assumption w(k) = w o = w^, 

since the frequency spectrum extends over a large band Aw > w above 
0 )^ for y/Z << 1. For these reasons, the free energy proposed by them 
is not applicable to proper nonideal plasmas, y/Z ' 1, nor is it correct 
for weakly nonideal conditions, y/Z << 1. 

2. In termedia t e Nonideal Plasmas. For intermediate nonideal 

conditions, 1 . y/Z < 10, the spectrum m(k) of electron oscillations 

extends over a region Aw 0[w ] above w 

P P 

kr 1 k r ~ 1. Also in this case, a relatively simple formula can be 

e e t 

devised for the free energy. The logarithmic term in F , Eq. (25) 

is negligible compared to that in F , Eq. (38). for y/z'■ 1 since 

t .-.Tic iSk./KT for y/Z > 1 bv Eos. (45) and (46), respectively, 
p Mi 

Accordingly, the nonideal part (48) of the free energy is for inter¬ 
mediate non ideal plasmas: 

E/NKT = -7ty + (3/Z)fuy + (3/Z) f« (8.^/v^ - (3/Z)0(r 
- 3f : (t ,ak (j ) , y/z T 1. 


by Eq. (6) 


(n/n) 


( 53 ) 
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For )/'/. - 1, the ions can bo assumed to bo non-do^cno rato , t ^ - fic^ok /KT ■ 1 


by Eq. (4b), so that the ion integral (41) rodnoos to 

«(>'.) = i/i. r. 1 . 

Since t ^ r 1 and ak^t r . 1 [ Kq. (A r ,)] for y/Z ' 1 > ti, t . electron 
integral (U)) is significantly smaller than l ; (r .) = 1/3, 

0 • I (, ,ak )• (z 2 ~z ? )^ 2 (ak e ) ~ J f f f (I-c*" e °/1-e P °) -- 1, y/Y.i 1. 
peep e p 1 

The lower and upper bounds of F(i_^,ak ) have boon obtained by means of 

the mean value theorem for the integral (30), 


((c ,ak ) = (ak e )' 
p e e p 


V ^ 


-l)d> 


(14) 


(05) 


( 56 ) 


While for strongly nonideal conditions, the contribution of the 

electron oscillations to the free energy is completely negligible, 

this contribution is still insignificant for intermediate nonideal 

conditions, y/Z > 1, by Eq. (55). For more exact evaluations, the 

small term F(c ,ak ) in Eq. (53) can be computed from Eq. (30) or (32). 
p e 


3. Weakley Non ideal Plasma s. Although the theory of weakly 

1-5) 

non ideal systems is well understood, it is interesting to investigate 
wheiher the present model for proper nonideal plasmas gives reasonable 
results in the limit y/Z << 1. For y/Z << 1 it is ak^ >> 1 by Eq. (45), 
and the spectrum w(k) of electron oscillations extends over a large 
region Am '•> u> above w by Eq. (5). The electron integral becomes for 


ak > > 1, 


F(e ,ak ) = c 

P e P 


a k 

(ak e )" 3 J (e C P X -l) _ 1 x 3 dx, y/Z « 


1, 


(57) 




r( , ak o )='/ i | , -l( 1 ak^) F i(r ak p ) 2 - ...], r .* « 1. (58) 
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.Aititougli t ak^ is independent of y/Z by Eqs. (27) and (23), the expansion 
(58) is valid since the electrons are certainly nondegenerate, }■ / r <"•' ] 


for y / 7. << 1, and 

1 / 1 / 1 / 

e ak = ( 4 ) • (9n 2 /4) 3 X /r 1, > / r <• 1. 

p e e e e e e 


(59) 


For nondegenerate ions, the integral (41) is C(t^) = 1/3 by Eq. (43) 
since e. --v 1. Thus, one obtains from Eqs. (18), (25) and (38) for 
the interaction part of the free energy ol weak 1v nonideal plasmas: 


AF/NKT = -u(y)y + (3/Z)i';ty + (3//.) tn (dc /v ) 

+ 36l(e p ak e ) - (1+Z _1 ), y/7." 1, 

where t.bo logarithmic term in Eq. (25) has been o- ponded for 


(60) 


c ak ■* 1. 

P e 

In Eq. (60), u(y) is the Madelnng constant of the weakly nonideal 
plasma with weak electron and ion ordering, n(y) *0 for y -» 0. Comparison 
of the term -a(y)y(NKT) in Eq. (60) with AF = 

- (NET) (2/3) it ^ (1+Z) ^e 3 n ' ‘ (KT) ^ of the Debve-Hueckel theory—^ 

(weakly nonideal plasmas) yields the result 

<i(y) = (2/3)ir !/> (l+Z~V /? y l/2 , y/7. < 1 . (61) 


The previous theories of weakly nonidoal plasmas do not lead to the 
logarithmic terms in Kq. (60) since they ’o not take into account the 

effects of electron and ion oscillations. 

The presented theory is applicable to nonideal plasmas in the gaseous 
phase (0-y-l) and the quasi-liquid (l-v-y,) phase. Whether the Eqs. (25) 
and (38) are applicable to hot (T>10’ C K) plasmas in the solid phase 
(y y ) cannot be judged at this time, since not enough is experimentally 
1 own about the latter, extreme state ol matter. 





APPENDIX: Expansion of I (t. ,ak ) 

1> <• 

The integral (23) is conveniently rewritten in the fori: 

f(e ,ak ) = (ak e ) I (e ,c) 
p e e p p 


where 


l V‘> j < • > -e) % (ee 1 )- ; J< , o- ,. p 


Since r > 0, i.e. e < 1, there exists the series expansion, 

(e*-1) 1 = £ e mC , e > 0 . 

m=l 


The snhst itntion, u =c - e . du = dt , and Kq. (A3) tiansiorm Kq. (A2) to 


Kr,,r) = l e m P 


me r V h -mu 


r •/ 

I u 


(u+2e p ) (• du. 


For u< 2e , i.e., e < 3e , the binomial expansion, 

P P 

(u+2e p )^ - (2c p )^ l u/2c p <l, (/ 

n=0 

is used, which reduces Eq. (A4) to the double series: 

I (t- ,e) = (2 l ) l e“ m 'p l ( n ) (2r )~ n m _( 7 +n) Y(4-+ n, (c-e )m) , 


m=l n~0 


c < 3r , 

P 


where 


5 . e-e 3 . 

-z- + n . p -rr+n 

5 , , , 2 f 2 - |nu , 

V (—+n , (e -c )m) = m i u e du, 


is the incomplete gamma function, which is tabulated. In an 


amologous way, the integral (A2) can be solved for u > 2e p , i.e., 
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V. Variational principles and canonical theory 

FOR MANY - COMPONENT PLASMAS WITH SELF-CONSISTENT FIELD 


ABSTRACT 

The Lagrangian L «* /// L (q , Dq ,/Dt, Y*q )d 3 r and Hamiltonian 

8 S 8 

H /// N(p , q , V-q )d 3 r for many-component plasmas with self-consistent field 
s s s 

interactions are derived. The canonical fields p (r,t) and q (r,t) are defined 

S 3 

by p = m v and Sq /St = n v , respectively, where n (r,t) is the density and 

S 8 S 3 S 8 S 

■ ► —► 

v (r,t) is the velocity of the component s, s * 1,2,...N. Based on the action 
principle, the Lagrange and Hamilton equations of motion for the components (s) 
are presented as functional derivatives of L and H with respect to the canonical 
momenta p and q^. it is shown that the new formulations oi many-component 
plasma dynamics have mathematical advantages compared with the conventional 
variational principles and hydrodynamic equations. 
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INTRODUCTION 

Since the work of Clebsch — it Is known that the space and time dependent 
dynamics of perfect fluids can be derived from variational principles. - 
To-date, iiowever, no variational principle has been found the extremals of 
which give the Navier-Stokes equations with viscous dissipation. A comprehensive 
review of the most significant contributions to variational principles in fluid 
mechanics has been presented by Serrin. — Similar variational principi s for 

perfect one-fluid plasmas and magnefohydrodynamic fluids have been given by 

7 8 

Newcomb — and Zakharov. — 

Herein, a new variational principle and canonical theory for perfect, non- 

relativistic many-component plasmas is developed, in which the charged particles 

interact through the self-consistent electric field. Interactions through 

.••elf-consistent magnetic fields are neglected for non-relativistic velocities 

v i. and temperatures kT << m c 2 of the electrons. We show that the electro- 
e e e 

hydrodynamic equations of the components s = L, 2, ... N of the plasma are 

9 

derivable as functional derivatives — of the Lagrange function 

q -V , --V 4 1 

L = /// L(r,t)d J r and Hamilton function H = Iff n(r,t)d r of the plasma of 

it 

volume il. The Lagrange L and Hamilton tt densities are functionals of the 

"► ">■ . 

canonical congugate variables and defined by = m^v^(r,t) and 

8q /Dt = n (r,t)v (r,t), where m is the mass, n is the density, and v is the 

S S S S S s 

(local) average velocity of the particles of the s-th component. 

The variational principle and canonical theory for many-component plasmas 
provides a theoretical foundatIon for the analysis of classical phenomena such 
as nonlinear electron and ion acoustic waves, space charge fluctuations, 
electrostatic turbulence, and electrohydrodyn.imic initial-boundary-vn1ue 




problems. The Hamilton function and canonical formalism were derived as a 
theoretical basis for a statistical mechanics of hydrodynamic plasma fields 
concerned with the evaluation of the free energy of electron-ion plasmas, which 
are nonideal due to the sel f-consist ont field interactions. The statistics of 
fields and wave modes represents an alternative to the treatment of many-body 
systems within the framework of the discrete particle picture. 
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VLASOV PLASMA 


The subject of consideration is a plasma consisting of an electron component 

and N-l ion components, which interact through the self-consistent electric field 

E(r,t) of all charged particles. The hydrodynamic equations for each component 

s of the plasma are obtained as moments of the Vlasov equation for its velocity 

distribution function f (v,r,t), and are in absence of collisions (s=l ,2 ,. .. N) 

s 

7 — (n m v ) + V■ (n m v v ) = - VP + n e E, (1) 

3t s s s ssss s s s ’ 

|— n + V-(n v ) = 0 , (2) 

31 s s s 

|- P + V• (P v ) = -(-> - 1 )P V-v , (3) 

Sts ss s ss 


where 


V-E = e " 1 V n e (4) 

o u , s s 
s=l 

-> ■> ->■ 

relates the self-consistent field E(r,t) to its space charge sources n (r,t)e . 

s s 

(c is the dielectric permittivity). The charge of the particles is e = -e for the 
o s 

electrons and = Z^e, |z | = 1 , 2 , 3 ,..., for the ions (e is the elementary charge), 
The moments are the density, velocity, and pressure fields cf the s-th component. 


n g (r,t) = /// f g (v,r,t)d 3 v , 

—<X) 

-foo 

v (r,t) = n 1 /// vf (v,r,t)d 3 v , 

s s s 

—00 

+°° 

P (r,t) = (2/3) III %m (v -v ) 2 f (v,r,t) d 3 v. 


P, v ... 


Furthermore, m is the mass of the s particles and = c /c is the 

S & So 

adiabatic coefficient of the s component (v is the particle velocity). 

Equations (1) - (4) describe the macroscopic dynamics of an N-component 
Vlasov plasma, — in which external magnetic fields are assumed to be absent, 
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and internal magnetic fields and binary correlations are negligible. 
Maxwell's equations indicate that in the absence of magnetic fields, H = (5: 


V x E = 0 


-r 

;»E v -* A 

— + ) n e V =0 

It . S S S 


Thus, the self-consistent field K(r,t) is longitudinal, and the displacement 

(e 3E/3t) and convection (n e v ) currents annihilate themselves at every point 
o s s s 

(r,t). This effect was first noted by Ehrenfest in his discussion of 
"accelerated charge motions without electromagnetic radiation". — Typical 
examples for accelerated charge motions without transverse radiation field 
are electron and ion acoustic waves, and space charge turbulence. 

Substitution of Eq. (b) into Eq. (4) yields the 
conservation equation, 0 £ ^e^/lt -f V-^ne^v = 0, fcr the space charge. The 
latter equation is, therefore, not independent [in addition, it is obtainable 
by summing Eq. (2) multiplied by with respect to s). 

Due to the absence of i) external magnetic fields and ii)binnry interactions 
(vanishing viscous and thermal dissipation ), the pressure in Fqs. (1) and (3) is 

isotropic. The closure of the Eqs. (l)-(4) is accomplished by setting all 

higher order moments of the velocity distribution zero, such as the off- 

diagonal stress components, n„ ~ 0, the third order heal flux tensor, 

q. - 0, etc. The resulting electrohydrodynamic or hydrodynamic Vlasov 
t J k 

equations are generally accepted in the literature. 
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EXTREMUM PRINCIPLE 


The dynamics of the many-component plasma is now condensed into a varia¬ 
tional principle for the Lagrange integral functional L = /// L(r,t)d 3 r, the 
extremals of which give the Eulerian equations (1) for the components s = 1,2,...N. 

In this approach, Eqs. (2) - (A) are auxiliary equations by means of which the 

^ 

fields n g (r,t), P (r,t), and E(r,t) can be eliminated from Eq. (1). This 
elimination will be carried through after introducing canonical variables. 

In analogy to the Lagrange density L(r,t) for nonconducting fluids,— the 
Lagrange density for the many-component plasma with the self-consistent field 
E(r,t) is sought in the form 


L( r,t) 



_> 2 
n m v *■ 
s s s 


Y 




(7) 


If Eq. (7) is the correct Lagrange density, the Eulerian equations of motion 

6 

for the plasma components should follow from the extremum condition — 


<S/dt///f(r,t)d 3 r = 0, 


( 8 ) 


r. e., 
t„ 


N 


fdtfff ) [n m v -dv + — m v ? 6n --- 

i j j ‘ s s s s 2 s s s v - 1 
t 2 s—1 's 


6P 


- e E* 6E]d 3 r = 0. 
o 


(9) 


The Eulerian virtual displacement, O y (r,t), of a mass element of the s-th com¬ 
ponent, 


► 

o 

s 


(r,L) 


v 

s 


(r,t)6t, 


( 10 ) 


vanishes at the initial (t^) and final (1 2 ) instants of the .variation of the 
trajectory of the fluid element which occurs in the interval t^ < t < t^» 


a (r,t ) = 5, a (r,t ) = 5, 

s 1 s £ 


( 11 ) 
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and is tangential to tne planes dA bounding the plasma at r = s: 


0 s (r,t) -dA(s) = 0. (12) 

The variations of the fields of each component s are given as functionals of 
its displacement o^Cr.t): 


* 

iSv - 

3 o -»■■•>■ -»■ 

s + v "S/a - a *Vv , 

. ^ s s s s 

(13) 

s 

31 


= 

- V. (n o ) , 

(14) 

s 

s s 

6 P = 

- o -VP - y P V*o , 

(13) 

s 

s sss 

~y 

f N 

— 1 n -V 


<5E = 

-tinea 
o ^, sss 

(16) 


Equation (13) is obtained by setting the variation of the Lagrangian velocity 
equal to the substantial derivative of the Eulerian displacement. Equations 
(13)—(16) result directly from Eqs. (2), (3),and ( 6 ), by means of Eq. (10). 
Substitution of Eqs. (14)—(15) into Eq. (9) gives 




-> 

Do 


/dt/jf ) [n m v • (-r— + v • Vo - o -Vv ) - ~r m v z V-(n o ) 
1 J J L , s s s 91 s s s s 2 s s s s 


t U s= 1 


+ ——:—(o *VP + y P V*o ) + n e E-o ]d 3 r = 0 
Y - 1 s s ss s ss s 


(17) 


By means of partial integrations. Gauss' integral theorem, and the surface 
condition (12), Eq. (17) is reduced to 


t.. 


N 


fdtjjf ) a • [— (n m v ) + V*(n mvv)+VP - ne Eld 3 r = 0 . (18) 

. K 3t sss ssss s ss 

t SI S= 1 


Since Eq. (18) is zero for arbitrary variations a , s = 1,2,...N, the integrand 
in the bracket must vanish a‘: any point (r,t) of the plasma. Thus, the Eulerian 


equations of motion (l) are obtained for the N components from the variational 




principle (8). This completes the proof that liq. (7) is the correct Lagrange 

density for the multicomponent plasma. 

The Lagrange density L and the Hamilton density If are not only functionals 

of fields of r and t, but also depend, in general, explicitly on r and t. In the 

-y 

following canonical theory, the explicit r,t dependence need not to be formally 
indicated, since the virtual displacements constructing the varied paths are for 


constant r and t. 
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LAGRANGE EQUATIONS 

!n the Lagrangian formulation of many-component plasma dynamics, the 

Lagrangian L =/// L(r,t)d 3 r is expressed as a functional of the fields q^Cr.t), 

-y ~y 

which will be shown to form with the fields p g (r,t) canonical variables. These 
are defined by, s = 1,2,...N: 

p = m v , 3q /3t = n v . (19) 

s s s s s s 

■+ -> -y 

Note that p g is proportional to the velocity field V g , whereas q g is a t-integral 
functional of the particle flux n s v s > Integration of the auxilary equations (12), 
(3), (6) (Eq. (6) is equivalent to Eq. (4)], and Eq. (19) yield 


n = v - V *q , v = v (r) 
s s s s s 


P = P [(v - V*q )/n ] 

s so s s so 


> -> -1 -*• 

E = e - e e q , e = e (r), 
S OSS s s 


( 20 ) 

( 21 ) 

( 22 ) 


v = (3q /3t)/(v - V-q ) 

s s s s 


(23) 


P , n , v ,and e are constants with respect to the integration variables. 

SO so s s 

Substitution of Eqs. (20) - (23) into Eq. (7) gives the Lagrange density of the 


plasma as a functional of the q g , 
N ^ ^ ^ 

L = I L S ( V 3 q s /3t ’ v ' q = ) 


S=1 


(24) 


where 


i = % m 
s s 


(3q /3t) 2 P 


v - V*q 3 

s ^s 


v - V • q v c , 

so , s s . s o -I t v? 

-_ (-)-— ( E _ c e q ) / 

> - 1 n 2 soss 


(25) 


so 


is the Lagrange density of the s component. The Lagrangians of the plasma and 


i 

'i 
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its components are given by 


•L = I L , L = J/J id 3 ?. 


(26) 


s=l 


The fields p g and defined in Eq. (19) are canonical variables since by 


Eqs. (24) - (25) 
> 3L 


3q /3t 
s 


S 3(3q g /3t) v g - V-q g 


— = m v 
s s 


(27) 


Thus, by introducing the canonical variable q g , JL has been decomposed into 

the L g of the individual components. Note that the original Lagrange density in 

Eq. (7) does not permit such a decomposition due to the coupling of the components 

by the common, nonlinear field energy e F.^/2. 

o 

For the Lagrange density of the form of F.q. (24), the variational principle 


in Eq. (8), 

r 


N C 2 


/ dt Jff Ld 3 r = l 6 / dt /// L d 3 r = 0, 


s=l t„ 


(28) 


becomes, since the L g are independent for s ~ l,2,...N, 
1 2 

<5 / dt [ff L (q , 3q /3t, V’q )d 3 r = 0. 

* is s s s 

t J u 


Hence 

t.. 


3q 3 L 


f dt Jff l~-5q + --•6( 3 - t -%+ —V 6V-q ]d 3 r * 0 

t, n 3q 3(3q /3t 3 V* q 

1 s s s 


(29) 


(30) 


where 


6 q = (3q /3t)6t=nv6t=no , 
s s s s s s 


(31) 


<$q (r,t.) = 5, 5q (r,t ) = 5, 5q (r,L) -dA(s) = 0, 
si s 2 s r=s 


(32) 
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by Eqs. (10) - (12). Ry means of the partial integrations, 

31 3q 31 t., C 2 „ 31 


l 2 31 3q , 31 t C 2 

/--A- -6 (^-)dt = [-—--6q 1 

t 3(3q s /3t) D(3q s /Dt) t, t 


- / 3q_- ±-dt, 

3(3q g /3t) 


1 


(33) 


ill 

S2 3V*q 


3L >1 

—3- 6 V • q d r r = i -2- ;q 

D V ■ c] 


- /// aq 0 ‘ v ( 


31 


Q 


s 3 V • q 


L -)d 3 r, 


(34) 


and Eq. (32), the extremum principle is reduced to 
1 2 3l , 31 31 

/dtj//l-^ - £ - v(~A.)].^ d 3^ = o 

t, 3q 3 ( a q 3t) 3V-q 

1 *s S S 


(35) 


The condition for the vanishing of this integral for arbitrary 6q g leads to the 
Lngrnngian equation of motion of the s component of the plasma: 


3__ 
3 t 


31 


31 


31 


-T + V(—^-) = , 


(36) 


3(3q g /Dt) 3q g 37-q^ 


whore 


3q /3t 


w = m (- 
s s 


) X V x (- 


3q /3t 


(37) 


V - V'q 

s s 


v - V* q 
s 's 


is a vector field perpendicular to 6q , w ■ 5q = 0 by Eq. (31). It is noted that 
the extremals of Eq. (35) are undetermined up to fields 

> -> -* > •r . - > ■> . 

w i iSii (w . m v x V y v 7 0 for v f -Vi p ). 

s ssss s s s 


In terms of the I.agrangian of the multi-component plasma, Eq. (26), and 
9 

Junctional derivatives,— the Lagrange equation (36) assumes the most general 
form, 

3 jI, 61. 


1 6(3q /3t 5q 

s £ 




(38) 


where 


41, 

-> 

5q.. 


31 

-r- 7 (—-J- 

3,, a 


6L 


31 


3 V • q 


6(3q /3t) 3(3q /3t) 

S S 


(39) 
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Subst i t ut ion of L f tom Fq. (25) into Kq. (36) indeed yields the fc 


iquation (l). 


.s into V (v 1 ) 
s 


^ - + 

v * Vv 
s s 


i 


1e r Ian 
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HAMILTON EQUATIONS 


Another important formulation ot the dynamics of a many-component plasma is 


obtained by means of the Hamilton forma 1 i ;;m. 
momenta, the Hamilton density of a component s is 
3q 

H = p ~ - L , 

s s Dt s 


Since p and q are canonical 
s s 


. p v V • q , c , „ 

H = -— (v - v-q ) P^ + - 7 (-) + -7 (t - e e q ) (41) 

s 2m s s s Y - 1 n 2 s oss 

s s so 

by Eq. (25). Accordingly, the Hamilton density of the plasma is a functional of 
the form 

N 

H = 'i H (p , q , V-q ) . (42) 

u s s s s 

S=1 

The resulting Hamiltonians of the plasma and. Its components are 

H = [ H H = /// H d3r. (41) 

s—1 b s n 

Thus, also within the Hamilton approach, a decomposition of the Hamiltonian 

H into component Hamiltonians H is achieved, each 11 depending only on the 

s s 

fields of the component s. This is quite remarkable, since the various components 

r > 

s Interact through the self-consistent field E(r,t). The decomposition is 
physically possible, since the electric fields of the different components 
superimpose in a linear way. 

The Hamilton equations of motion of the components s follow from the action 
principle (29), 


<5 / dt /// (P/T,:* - = °. 


(44) 







which gives, since the Hamilton density is a function of the form 


H s = H s (l v v V ‘i ) ’ 

*~2 Dq 3H ^ 3q a 1/ jH 

! dt IS! Uj t - —r>' 6 P a + P s ' 6 ( 3T ) -- i“ *'■’*%]d r = 0 . (45) 

t, ft 3p ‘ oq jV•q 

I s s n s 

Interchange of the sequence of variation and differentiation and partial inte¬ 


grations , 


2 , - > 3 


t., 

t 0 l . v - 

2 f ox 


/? 8 *F q s dtB iv\i t - ri± -*vit. 


l i ■ 1 4 :)l 

3tf "jH dH 

US - - V-<5q\d 3 r = if -~ <Sq *dA - /// 6q -V(--)d 3 r, (47) 

ft 3V.q. b 3V-q, S ft S 3V-q 

b S S 

where the variation 5q^ satisfies the condition in Eq, (32), transform Eq. (45) 


t 2 9q 3 H 3p 3 H 3 H 

1 dtj// ~ ~)’6p - + V(—S-)].«q >d*r = 0. (48; 

C 1 « 3 P S 3 ^ s 3 ^ s 

The conditions for the vanishing of this integral for arbitrary variations 

-> -V 

4p^ and 5q^ gives the Hamilton equations of motion for the plasma components s: 

3q 3H 

s s 

IT = ~ . (Ay; 

3p s 

3p 3 II ' 3 H 

s s „, s . -> 

*T +— - v(— —) = w g , (50) 

3q __ d7*q 


where 


-> 1 >■ -> 
w = — p x V x p 

*■’ TTI ft ft 


m s 
s 


is a vector field perpendicular to iq = n v ot, i.e.,w •6q =0. As above 

s s s s s ’ 

w s is determined from the condition of rotational flow (w f ft for V x v / S). 
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By introducing the Hamiltonian ni the p!asina , Eq. (ah), and functional 

Q 

derivatives,- the Hamilton equations (49) - (50) are brought into the most 
general form. 



at 


> 



(52) 



3t 


4H 


(53) 


where 


6H 

* 

OH 

s 

-► ’ 

V jH s 

Oil b 

51/ 

- V (— J) 

(54) 

" P s 

Dp 

s 

fiq Dq 

s s 

3V-q s 



Evaluation of Eqs. (49) - (50) for the Hamilton density in Eq. (41) results in 
Eq. (19) and the Eulerian equation of motion (1), respectively. 


CONCLUSIONS 


Comprehensive Lagrangian and Hamiltonian methods have been developed for 

many-component plasmas with internal interactions through the self-consistent 

electric field. The theory is based on the Lagrange density (26) and the 

Hamilton density (41), and the canonical variables p^ and q defined by 
' ► -> . '► 

p = in v and 3q /Dt = n v . respectively. In the limit of a single component 
b s s s s s 

(N=l) and vanishing Coulomb interactions, the corresponding Lagrange and 

Hamilton equations of motion for an ordinary neutral particle gas result. 

The canonical theory presented reduces the coupled partial differential 

equations (1) - (4) for the fields v (r,t), n (r,t), P (r,t), and E(r,t) to 

s s s 

a single partial differential equation for the canonical field q g (r,t). The 
fundamental Lagrange equation of motion is 
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mathematical advantages of the theory presented are obvious, in particular for 

plasmas with several components. In comparison with the previous formulations 

13 4 5 

of variational principles in fluid mechanics —*—*—*— which make use of an 


excessive number of scalar Clebsch parameters, the present canonical theory in 

-+ -> 

terms of the canonical vector fields p g and excells in mathematical simplicity. 





69 


Referen ces 

1. A. Clebsch, J. Mathematik _54, 38, (1875). 

2. H. Bateman, Proc. R. Soc. London Ser. A 125 , 598 (1929). 

3. L. W. Herivel, Proc. Cambridge Philos. Soc. 51^, 344 (1955). 

4. H. Ito, Progr. Theor. Pliys. 1_3, 453 (1955). 

5. F. P. Bretherton, J. Fluid Mech. 44^, 19 (1970). 

6 . J. Serrin, " Mathematical Principles of Classical Fluid M ec hani cs", 
Handbook of Physics (Springer, New York, 1959), Vol. ill, Pt. 1, p 

7. W. A. Newcomb, Nucl. Fusion, Sttppl., Pt. 2, 451, (1962). 

8 . V. E. Zakharov, Zh. Eksp. Teor. Fiz. 60, 1714 (1971), [Sov. Phys. 

JETP 33, 927 (1971)]. 

9. J. Rzewuski, Field Theory , (Hafnei, New York, 1969), Pt. 1. The 

functional derivative 6 F[¥(x)]/6Y(x) of a scalar functional F = 

F[f(x)] is defined by [ 6 (x - x')= Dirac 6 -function]: 

5F[ Y(x)] F[Y(x') + Y 6 (x - x')] - F[Y(x')] 

. v — = lim —-———--— — 

6nx) e ->0 

10. A. A. Vlasov, Many Particle Theory and its Applicat i on to P lasma 


in 

144. 


(Gordon and Breach, New York, 1961). 

11. P. Ehrenfest, Phys. Z. LI, 708, (1910). 



70 


VI. Distribution function of turbulent velocities 


ABSTRACT 

The Hamiltonian H and Hamilton equations of motion are derived for 
-► -► 

the Fourier amplitudes p£ and q+ of the canonical conjugate fields p and 
•+ -> ■+ “► -► 

q defined p = mv and Dq/3t * nv, vhere v(r,t) and n(r,t) are the velocity 
and density fields of ideal, compressible gases in the state of fully 
developed turbulence. A Liouville equation is presented for the 
distribution function f(p^,q^;{lc}) in the multidimensional phase space 
formed by the scalar components of the set {£} of wave mode vectors 
P£ and q^. As an application, that stationary solution f « f(It) of the 
Liouville equation is calculated, which maximizes the turbulence entropy. 

It is shown that the distribution of the velocity and density fluctuations 
of compressible gases is Gaussian in fully developed turbulence, in 
agreement with the experiments. 









INTRODUCTION 


The turbulence probLem has aroused the interest of many celebrated 

1 ) 2) i) 4) 

researchers, such as Reynolds , von KArman , Taylor'', Heisenberg , 

Kolmogorov , Prandtl^, and Frenkiel^, to name only a few. In spite 

r ,, , 8-10) , L 

of all progress made to date , many experimental observations on the 
turbulence phenomenon are still unexplained from first principles. One 

of the most elementary experimental results is the Gaussian velocity 

-►2 -2 -> 
distribution ~ exp(-v /c ) of the turbulent velocity fluctuation v in 

11 ) * 

compressible gases. The Gaussian velocity distribution of fully developed 

turbulence is a fundamental fact as the Maxwellian distribution 

1 ->2 ■> 

~ exp(—^mv /KT) of the molecular velocities v in statistical equilibrium 

gases. 

The Maxwellian velocity distribution, which is the result of many 
molecular collisions, is independent of the molecular interactions, which 
have brought about the statistical equilibrium. The experiments indicate 
that in fully developed turbulence, the distribution function of the 
velocity fluctuations shows only in the far-out wings jv ]-*• « smal 1 deviations 
from the Gaussian distribution.^^ Apparently, the distribution of the 
velocity fluctions in stationary turbulence depends hardly on the viscous 
momentum transfer, which is negligible for the velocity fluctuations of 
large scale but significant for the velocity fluctuations of small 

scale A-> 0 , concerning the sustainment of the turbulent state since 

-V^v,*- ~ (A = 2-n/k) . 

k k 

In other words, the Gaussian distribution of the velocity fluctuations 
in fully developed turbuLence is mainly sustained by nonlinear mode coupling 
of the turbulence elements of different wave vectors ic. For this reason. 


it should be possible to derive the turbulent velocity distribution from 
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a statistical hydrodynamics of the turbulent llnid notions, without taking 
viscosity into account. By Fourier analyzing the canonical conjugate 

- V - ► -+ 

variables p(r,t) and q(r,t) of the turbulent gas, it is shown that the 
-> -> 

Fourier amplitudes p^(t) and q£(t) satisfy Hamilton equations of motion, i.e. 
the random wave modes k behave like quasi-particles. This makes it possible 
to derive a Liouville equation in the hyper-space lormed by the scalar 

• V -r 

components of the set lp^ , q^} . From this fundamental equation, the 
quasi-Gaussian distribution for the turbulent velocity fluctuations is 
derived as that stationary solution which maximizes l he- entropy of 


turbulence.!. 




HAMILTON EQUATIONS 


In the statistical turbulence theory to he presented, tiie effects of 
viscosity are disregarded so that the random velocity field v(r,t) of the gas 

- t- —V 

Is irrotational, Vxv=0, in the absence of external forces which could 

> > 

artifically produce velocity curls. The condition Vxv=0 defines longitudinal 

or so-called "acoustic" turbulence in compressible gases.It represents 

a good approximation for the treatment of certain turbulence phenomenon, 

such as the present calculate on of the turbulent velocity distribution. 

The basis for the continuum-statistical considerations form the Hamilton 

equations of motion for the compressible, nonviscous gas. If v(r,t) 

• ► 

and n(r,t) designate the velocity and density field ol the gas, respectively, 

■v > ir ■ v 

the canonical conjugate fields p(r,t) and q(r,t) of the gas are defined by the 
proportionality (m=mass of gas molecule) and functional relations,respectively 


p = mv 


3q/3t = nv 


The Hamilton equations of motion and the Hamiltonian H of the gaseous system 

of volume ft = fff d *r = L L L are given by:*^ 

x y z 

3q _ Mi 

3jV = _ 6H 

3t 6cf , (3) 


, 1 V ^ () / v~ V * 

H = -M’q)p + —f — 


= 3H 

6 p 3 p 


fib = _3H _ „ /_QH_ > 
6 3 ~ 3d ' \3V-q 


delines the functional detivatives 4 |3H/3q 0 since If = H(p,V-q) by Eq.(A)|. 


In Kq.(A), y (, = (■/<•) is Lite aiMabati,- coefficient, P and n 

o o 1 V o o 
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[reference values of P(r,L) and n(r,t)] and v are constants obtained through 
integration of the polytropic energy conservation and continuity equations of the 
gas, respectively. 

Since n(r,t) and v(r,t) are random turbulent fields, p(r,t) and q(r,t) 
are represented as complex Fourier series, 

P (r, t) =£Pjj(t) exp (ik-r), q(r,t) =^q^(t) exp (ik-r), (6) 


where 

k = 2 tt 


V M \> I 

x _ z\ 


L ’ 1. L 
x y z 


{k} 

= 0, ± 1, + 2, 


, ± v 


(7) 


x,y ,2 

are the discrete eigenvalues (Q<°>) of the wave vector k. Since the minimum 
wave length X = 2n/k^ is of the order of the mean free path, Z, one has 
v ~[L/£]>> 1 for Z<< L. 

Hamiltonian equations of motion similar to those for the canonical- 

■ V V -k - F 

conjugate fields p(r,t) and q(r,t) can be derived for their Fourier amplitudes 

V —V 

p,>(t) and q, k (t), i.e., for the individual wave modes. By Eqs. (2)-(b), 

q,i/dt (8) 


311 


Dji 


-k 


t) , i.e., tor the individual wave modes. By Lqs. 

for 


!l 


52 


and 


4 ■ f te 


57 


3 H 


3? 



+v . Mi. = 


I- 


ik ' r V'[-Jl d 3 


- 52dp//dt 


(9) 


< ► 3 -v ij<. r 

since Jj/7-A d r = 0 if A is expandable in orthogonal functions c = unit 


tensor). 

The canonical equations for the canonical-conjugate amplitudes p^vCt) and 
qj>(t) of the individual wave modes k are by Eqs. (8)-(9): 
d P2 


1 3H 


dt 


51 3<t 


( 10 ) 


-k 


and 


dq. 


1 311 


dt J 


(ID 


-k 
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Equations (10) and (11) represent a complete system of coupled nonlinear 
differentia] equations, which are antisymmetric in the k indices. These 
canonical equations determine the temporal development p'(t) and q + (t) of 

K K. 

>■ - r 

-► i k * r 

the wave modes k of spatial structure = e 


In Eqs. (10)—(11) the canonical equations of the wave mode dynamics are 
stated in complex form for the conjugate variables 


k = 


k* k = real 

(12) 


-> -V 

-y 


Si 

1! 

Q k + 1S '£ * 

Q k % S- = real , 

(13) 

so that 

> 


-y . *v -y 


P k = 

2 (p k + p -k } * 

k = 0/2i) (p£_p_£) » 

(14) 


2 (q k + q -k ) * 

-y ~y - *■ 


^k = 

S £ = (l/2i) (q k -q_ k ) , 

(15) 

since 


-y 

—rh -y 

(16) 

P £ = 

p -k’ 

q k = 

for real 

vector fields 

p(r,t) and q(r,t), Eq.(6). Substitution 

of Eqs . (12) 

(13) into 

Eqs. (10)-(11) yields the Hamilton equations for the 

real and 

imaginary 

parts of the 

canonical conjugate mode amplitudes: 


dP,>- 

k 

1 3H 

dP k 1 Dll 

(17) 

dt 

" ' n ’ 

dt Q 3?- 

k 

d % 

1 3H 

d?5>- , ... 


dt 

k 

k _ 1 311 

dt SI 3l+ 

(18) 


li 





MOUVILEE EQUATION 


The existence of the canonical equations (10) and (11), which arc- 
analogous to the canonical equations of a many-hody system, indicates a 
quasiparticle behavior of the individual wave modes. For statistical purposes, 
let a large number (ensemble) of similar, turbulent gaseous systems be 
introduced (each containing a set {k 1 of wave modes or quas ipartic 1 es Ic) 
in the mul tidemensional phase space formed l.y the real, and imaginary parts 
of the Fourier vector components p^> . and q^ ^ (i = l,2, 1) of the entire set 
ikjdefined in Eq.(7). The number of gas systems, which have their phases in 
the volume element il ^ ^ (d p^d^q j>) at the point (p£, q^;{k}) of the phase 
space, is given in terms of the density in phase space, f=f(p>, qj>, t; {It}), by 

K K 


dw = t'(p£, q k r, t; > k I ) H (d pod q^) , 

{k} 


whe re 


‘ d V’ R i;• ■ d V%- (20) 

in accordance with Eqs. (12) and (13). The phase-space density satisfies the 
continuity equation in phase space, since the number cj f plasma systems in the 
ensemble is conserved. 


Ft + "(it) ' <v (k> £) * 0 • 


where 


.vfei ^ 

■ kl Z- l<it ' dt 


\ 

«> = IrP*' 


The canonical equations (LO) and (11), which determine the motion of any 


point in phase space, indicate that the phase-space fluid is incompressible. 
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;► 1 X—' 

7 {k}‘ ik} ■ ri ]C 

(ki 

because ol this property, the continuity <• (|u 1 1 inn (21) can he brought into the 
form of the so-called Liouvi-llo equation: 




This fundamental equation describes the temporal development of the distribution 
in phase space. ft extends statistical mechanics to gas continua . in 
random motion. 

The Hamiltonian U in Eq. (4) is an invariant for a nondissipative gas. 

The steady-state solutions of Kq. (2b) are functions of H: 

f = f 00 , (26) 


since 


21 

dt 


I V'/lL >->H _ah iLL \ af 


0, for f = f(II) 


(27) 


In applications, the compact complex firm (2b) of the hiouville equation 
is mathematically preferable. In terms ot real vector components, Kqs. (22)- 
(23) become 

bit ■ E 

{k> 


f^'j 

l d t 


dR,* 

_ k 

dt 


d ^k 

dt 


(IS; 


dt 


(28) 


( k) 




(29) 


Since the phase space fluid is incompressibIe,V^ j * 

suit st i t lit ion of Kqs . (28)-(29) into the continuity equation 

the hiouville equation in real notation: 

a i_ .»n_ b[ _ ri,i _ 

Vn k ' T *k ' ^ ^ 



= 0 by Eq. (24), 
(21) results in 


1 
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Tlic* solution of the Louville equation permits the determination of tli 
macroscopic properties of turbulent gases as ensemble averages. The 
statistical considerations are based on the quasi-partic]e character of th 
individual wave modes. 
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TllRBlII.KNCK MNTROPY 


The state (if stationary, fully developed turbulence is one of maximum 
probability W. Since the entropy of a turbulent system is a function 
S " fn W of its probability U', the entropy oi turbulence must, assume 
an extremum in stationary state. 

The entropy associated with the turbulent modes of the gas is defined 

-y -y y 

as the phase space integral over the distribution function f = f(p^,q-+;(k) ) 

K K. 


■H 


S = - K | ... |f fllf H (d 6 p->d fV q>) 

{£} k k 


('ll) 


where k is a dimensional constant. This definition of turbulence entropy is 
in complete analogy of the definition of turbulence for arbitrary nonoquili¬ 
brium systems. 

The idealized description of the gas as a continuum without collisional 
dissipation is reflected in the time - indepencence of the turbulence entropy 
By I’qs. (3 ), (21), and (24). 


dS 

dt 


cj + Inf) ^ + v 


<3t 




(kl 


v,-:-, ;f n ( 

lk} J Ik) 


, 6 - 


f ydVdV) = 0. (32). 
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VELOCITY DISTRIBUTION 

As an application of the canonical theory, the statistical distribution 

of the random velocities in turbulent, compressible gases is calculated. 

Let the turbulent pas be, on Lhe average, homogeneous (n =■ n ) and at rest 

o 

(v - 0) so that for the spatial mean values 

p = 5 , q + d . 03) 

The gas is assumed to be in contact with an appropriate other system which 
sustains stationary, fully developed turbulence. The canonical conjugate 


fluctuation fields can then be expanded in the stationary Fourier series: 

p(t) = 2_J P k *f txplik-r I, q(r) = 2_j q^exp [ ik • r ], (14) 

{k} f k > 

where 


10 ) 


P k ' = ^ j"j'j*p<r)ox|»[-iE-r]d ( r , 




q(r)exp[-ik-r]d r 


(33) 

(3b) 


are the complex, time-independent Fourier amplitudes of the real fields p(r) 
and q(r): 


•>* _ -* 

P C = P -k’ 


-V 

q k = q -k 


CD 


The orthogonality relations for elementary vector fields of the form given in 
the above Fourier expansions are 

JJJ (P ) k‘ (q) k ,d r = !JJ P k' q lt' oxpfl(k “ k ') ,r l d ' r = u -CD 

Hy Eq.(27), the phase space density of a conservative gas (U=E) in a 
stationary state in on arbitrary function of the Hamiltonian 


. q£, fk1 ) = f[II(p£, q+; Ik} 


09) 


a i ua t state, tiie entropy of turbulence assumes an extremum(state 
r . bv l.q . (31 ) 







HI 


S(P£,q k -; ik I 


) = /r< P k ‘,qs fk) ) e» i(P k -.^; no ) :i(d V V ‘k 

1 k i 


) . (40) 


Ai cord inn! y > llie phase-space d ist r iliut ion 1(11) ol .stationary, lullv deve 1 oped 
turhuU’ine, is determined by the variational equation 

J. "* <’ti O li (d |),-d' <!,•I 0 , (41) 

where 


- - k I. . . I ( 1 I t’li 1) s f 11 (d ^|i -d f (| • I ■> 

i k } 


J...Jilll(d ,, 1 >' t ;d ,, q ^) = 0 , 

f. . . f H M i: (d'V/^q, ) - 0 , 

3 3 \i\ k " 


t > 


f S) 


since tin 1 number of jrik systems and the tot il < ■ 11 • rp.y i r e i. < 11 i. u. * . • ‘ | - 

ii.u 11 i p I i *-.it i mi i • 1 K*|. (41 ’) by ( <- 1 1 and Uq . < 4 ’■) l> > !■ , ■ ' < ■ i • ! in < • l 
iimJi t *■ mii nod kayraa;', i an mu 11 ip I i> r.-, i v*• ■ 

<■ r I ■ i-. • 

J. . .J < t W t + I + . It) l _U! p^-d 1^1 II 




1 + * ■* I! - 0 


i s 

1 1 ' ' i 1 I ■ ! 1 V 

. l-'i 'in ! <| . 

4 ■ ) tin 

rial 

i .’.it i i'ii 



(’ll 


[-.H|- (d*! 

• k :• 

t> - 

k’ 1 ’ V 


in t 1 ■ t iu ni 


f (H > 


■ k ■■) I 

/f 7 . , 4 • fi > 

• • -J I - •4i(fv,<|;; Ik. <) } [ (d p-d q^-). 

k 1 


(4; i 


[he Hami 1 toil i an in i '.< j . (4 ) becomes alt* r substitut ion **1 'be expansion 

i n Kqs . ()’>) — ( !b ) 

II 




14 H) 


- ‘ ^k-q >exp| i 


n "TT, k 
<) \ K I 


- k • 11 


.1 r 


( 49 ) 



-M 


since v 
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= n^. An evident bionomial expansion and subsequent spatial integration 
gives, under consideration of Eq. (38), 


where 


II - i:[- —ar— + ~ E pep,* + 1-yl’ n 2 E (ik'd 

L '" 1 2m un k k 2 °° :Ki 


k )( , k-q^ ) 


f h(p k ,q k’ <k ^ 50) 


.-1 



c-o' 


-J W n .'(1) .*(/-]), 

t U l 


E (ik n) - q y°i...ik v ’ f i , 1 j; it 


:('-i)■' (/-]) 

'k 


|k U) + ...k ( " 1) l-q_ (k -(l) + _ _ . .-(''-!) 


! 


■(k v " + ... + k' ij 
contains the terms ot third and higher order in tin Fourier amplitudes. The 

various sums in Eq. (3D) have the physical • -• *., i» i np. el spatial averages 

(t 1 . a e v. i I 1 s t lien I a in) 


151) 


V 


Vk = P 1 ,:| v 


'52) 


( ik ’<)/) ( i i "l|” 1 E'V " k * 5 

ik 1 


i > 11 


win i »■ v (v I') and u - n - u ate the : > h ■ 1 1 < • t • t veioc it',' util tlete.il'. 

t I it, I i la t i or. s a t the comp less i h 1 e re.. 

<‘e"ih i tiitij' el la|s. (\i) and (Mi) vield- the .t it inuarv i . t r i b u t i on 1 ii:i< t im: 

in tin- to mu : 


t (I! ) 


■ to mu : r __„ 

I ; ’oV • > * 1 - ’V • • 

• U..7fk' P k f 2' : ’o n o ^ Mk-q k 


I ( i k - q,* ) 


-.<h< P yq k *:<k*})J 


. >'* ) 


I e >••!>(;•■ I' (<-I) 1 I J. . expl-l ll(p^ ,q,-; k )] il d h p^d f *q^ 

i k i 


!,/ introducing the spatial averages defined in laps. (52)-(3 3) a physically 
more i I 1 ust rat i ve iopresfnt.it inn <>f the stationary turbulence distribution 
i s o!i t a i tied : 


I (II) 


A expj- r'..jj,n (| m v 2 + yyl’ 0 (n/n (1 ) 2 J - Fh(p^.q^; (k t) 


f 


(55) 






wile’re 


] = [oxpU-^JMy-I) 1 ) 1J*.. .Jexp[-3H(p *,qg; (k }) ) I! d 6 pod f 'q^ 


The parameter 3 of the distribution, which has the dimension of a reciprocal 
energy, is given in terms of the total r i i i l>n 1 < n t energy K contained in 
the I I netuntinns through Lhe relation 

p 

r /...j* IKp^c^tlf')! (plKpgoJ-tlkl) U d <> |» | _-ii , ’« I ‘ k - - ■ ( , ( ,1 ) 

Ik 1 

in v, i ew c>f the noi m a 1 i xat i on delined in Iq. !\h). f / < , — I j represents the 

o 

non t u rhu 1 eii L or thet.; il energy of tin- gaseous system of >. o 1 uii’e 

la|.(Si) indicates that t he d i s l r ilnit i on I n.i. t ,o,i f (M ) lor the turhulent 
1 1 in t u.it ions lias the tallowing prnpertii 


it Tlie d i s t r 11 ut i on . > 1 t he v > ■ I o, i t , ‘ I s. t u. ■ t i oie i 


" i an , 


. , nir • ( jH) 

I (v ) c 1 o 

i i ) Tin’ d i s t r i bu t i on o! the density : I u. t ua t i oils is o. cu.. . i an, 

. ,' 2. -\,V n . .«*’ . (V)> 

t (n > c- 2 o o 

ihe h i ghc r ctchr t c mis in the I i s t i i l.ut i mi I nm t I on i f 11 ) - • I I'q. (Vi ; I ue 


to h ( p • ,c| •, ; k •) , la |. ('» I ) . a t e ditti nit l 

e K 

lot low i tig reasons . I In v represent < i 


i vi' in i : i r intent •. I m tli 


• ■ I s o 1 eh -ns ; t s aid a- 1 :s■ i t 


t 1 ii> tii it ic is, i .. . they are uiiohsc i vab ’ • • it tin , /pet ' nr a 1 he ■ a ' >r r • its! 
only velocity or dens i l v I !m t it it ions. 'In. i , i m h (t ^ , c|, -; . I ) i.. snail el 
third and higher odor in the f 1 in t u.it i or amplitude; and is, in general, 
significantly smaller than the leading second order terms in the expansion 
"t H<p k >,>i k -; 'k 1) . 

Normal distributions of the forms given in 1‘cjs. ( r iH )-('>')) have been 
observed in fully developed l m i>; 1 ence of compressible gases.'Kr.per imental 1 y 
the distribution ol turbulent velocities is obtained by plotting, the Irequency 
of the occurence (over a larger time period) of a velocity signal with 
amplitudes between |v| and |v| + j Av| vs. |v|, i.e. by a time-averaging process. 



Ill the distribution function of Eq. (55), the quantities v 


and n 

represent spatial averages (Eqs. (52) and (51). This formal difference 
is, of course, no contradiction between experiment and theory, since the 
time average is equal to the spatial (ensemble) average in fully 
developed, homogeneous turbulence. In Fig. 1, the experimental (o) and theo¬ 
retical (-) distributions of the turbulent velocities ire compared!'^ 

It is seen that the theoretical (emission distribution of the turbulent 
velocities ! s in excellent agreement with tie e.-.per imental dale. 

The theory presented is highly idealized since viscous and thermal 
dissipation in t iie turbulent gas .ire not taken, explicitly, into account . 
it represents, however, a first attempt at extending statistical mechanics 
to r.indori continuous media, such as turbulent gases. Extensions of the 
theory to include viscous and thermal dissipation represent mainly 
mat hemal ical problems. 


i 
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VII. Appendix: stress relaxation naves in fluids 


ABSTRACT 

The Navler-Stokes equations for incompressible and compressible fluids 
are generalized by inclusion of viscous stress relaxation, as required by 
kinetic theory. Two initial-boundary-value problems of the nonlinear, generalized 
Navier-Stokes equations are solved analytically, which describe the propagation 
of transverse or shear waves due to temporal and spatial velocity pulses v(0,t) 
and v(x,0), respectively. It is shown that transverse perturbations propagate 
in form of a discontinuous wave with a finite wave speed due to viscous stress 
relaxation, whereas the conventional Navler-Stokes equations result in non¬ 
physical solutions suggesting a diffusion process covering the entire fluid 
with mlinile speed. 


! 





INTRODUCTION 


The nonlinear incompressible and compressible Navier-Slokes equations 
represent (quasi) parabolic and hyperbolic, partial differential equations, 
respectively. The former propagate signals with infinite speed and the latter 
propagate certain signals with finite speed in fluids. In an infinite, homo¬ 
geneous fluid, consider a small (linear) velocity perturbation, which is re¬ 
presentable as the Fourier integral, 

-f- -> -* 

'*/ f ►/IN iwt - ik*r , t 

v(r,t) - J v(k) e aw * 

— oo 

over elementary waves of wave length > = 2 tt/ k and frequency m(k) . If the 
fluid is compressible so that it sustains both pressure (p) and density (p) 

~ r, ~ 

perturbations, Dp/lt = c^ 3p/3t, the perturbation can propagate, e.g., 
in form of longitudinal sound waves with finite speed = (yP q /p ) 2 and 
dispersion law 

„? = c z k* + i(5>'/3p )mk ? 
s o 

In a fluid with a viscosity a perturbation may also propagate in form 
of a tiansverse or shear wave. If one applies the curl operation to the 
im 1 'inpi ess i hie or compressible Naviet -F>l"kcs equation, a dispersion law is 
to Mid tor the transverse perturbations whieb does not repieseut a wave 
phenomenon hut an aperiodie damping process witii 


= - (./p )k- 
o 


As is known, the acoustic dispersion law is derived from a hyperbolic wave 
eipiation, whereas the damping constant lot tho transverse modes follows 
from the parabolic vorticity equation (which is the same for incompressible 
and compressible fluids). From experiments, however, it is established that 
transverse perturbations (V x v^ = ik x v^ i tf) propagate as (hyperbolic) 


HQ 


shear waves with finite speed. 

It will be demonstrated herein that the (incompressible or compressible) 
Navier-Stok.es equations do not provide a correct description of shear waves. 

The discrepancy between the Navier-Stokes equations and the experiments on 
shear waves is resolved by introducing viscous stress relaxation, which leads 
to a hyperbolic transport equation for shear waves in incompressible or com¬ 
pressible fluids. For this reason, the transverse or shear waves represent 
"stress relaxation waves". 

As an illustration, two hyperbolic initial-boundary-value problems for 

shear waves with stress relaxation are solved. The solutions of the generalized 

Navier-Stokes equations with stress relaxation represent transverse waves 

which are discontinuous at the wave front and have a finite wave speed, 

c = (p/p t)' 4 < <« (t is the stress relaxation time). The first treats the 
o 

propagation of a shear wave into a semi-infinite fluid space, x - 0, produced 
by a temporal velocity impulse at the boundary x - 0 (accelerated wall). The 
second is concerned with the propagation of n shear wave into an infinite 
fluid space, - m - x - + caused by a spatial velocity pulse in the plane 
x = 0 at time t = 0. Both solutions are valid for nonlinear shear waves of 
arbitrary intensity. 


i 
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PHYSICAL PRINCIPLES 

In conventional fluid mechanics, —^ it is assumed that inhomogeneities 
V v in the velocity components v. produce instantaneously viscous stresses 

i- .1 i 

1! . Mathematically, this is expressed through the phenomenological "flux" ~ 

i j 

"force" relation 

II.. = - ,(V. v. + V.v. - ~ V v 6. .) 
ij i j j x 3 k k i.| 

f- 

where !• is the viscosity and 6 is the unit tensor. In a real continuum, 
however, velocity inhomogeneities do not switch on viscous stresses instan¬ 
taneously but in accordance with a relaxation process of characteristic time 
t. By means of the kinetic theory of gases -\nnd liquids, --^one can show 
that the transport equation for the viscous stresses has the form of a 
temporal (3/Dt) and spatial (v-V> relaxation equation. 


- - II. . + v, V, II . . 
Jt J.J k k jj 


t Si..- I i '(V.v. + V.v. - ' V v 'S . .) . 
ij l j j j i k k ij 


This equation is approximate insofar as the the coupiing oi heat flows <|^ 

and stresses II.. and higher order terms in the derivatives of v. are 
ij i 

2 3 ) 

neglected.—It has temporal and spatial derivatives as required for a 
J-L dependent field equation and is Cali lei covariant. If relaxation effects 
are disregarded, it reduces to the static stress equation. 

Thus, consideration ol viscous stress relaxation lends to a refornmlatI >n 
of the conventional Navfer-Stokes theory of incompressible and compressible 
Iluids. In place of the Navier-Stokes equations, we have the hydrodynamic 
equations with viscous stress relaxation: 

p( ft + ''• VV ' ) ~ V >' ' • (1) 


— + v.Vfi = - pV-v 

lit 


an > . -> 4-y , ‘ii 


+ v-vn + 


- = - - (Vv + Vv - 
x x ' 






Eqs. (1) - (3) hold for incompressible (V*v - 0) and compressible ( V • v j- 0) 

fit . For nonisothermal systems, the transport equations for thermal energy 

A 51 

and heat flux have to be added to Eqs. (1) - (3). —— 

If p and t can be treated as r-independent, It is mathematically more 
convenient to use instead of the tensor equation (3) the vector equation, 

— v-n + v* (v-vn) + t v-n = ~y.i (v 2 v + ~ w* v ), (A) 

since Eq. (1) contains the force density V-TT. If temporal and spatial re¬ 
laxation of the viscous stresses is disregarded, Eqs. (1) and (A) combine to 
the classical Navier-Stokes equation, p(3v/St + v*Vv) = - 7p + pV 2 v + (lj./3)VV* 
Equations (1) - (3) represent a hyperbolic system both in the compressible 
and incompressible cases. On the other hand, the conventional incompressible 
Navier-Stokes equations are parabolic. The corresponding field equations for 
incompressible fluids are obtained by setting V-v : 0 in Eqs. (2), (3), and (A) 


t > 
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INITIAL-BOUNDARY-VALUE PROBLEM FOR v(0,I)-PULSE 


A simple method for the generation of transverse waves in a viscous fluid 

consists in setting the plane x = 0 bounding a semi-infinite fluid (x _ 0, 

| s' 1 v >>■, I z I 1 "•) into sudden motion v = v H(t)e , where H(t) is the 
1,11 w y 

Heaviside step function. The resulting viscous interaction between the fluid 

and the accelerated wall produces a curl n x [y] - e v(x = Q,t) at the fluid 

surface which propagates in form of a transverse wave through the fluid in 

the x-direction. In this dynamic process, the fluid velocity is of the form 

v = 10, v(x,t), 0} so that V•v = 3v/3y = 0 and v*Vv = 0, i.e., the fluid 

motion behaves incompressible (even if the fluid is compressible) and linear. 

-> ■< —► —>■ ^—> 

Furthermore, vVJI = v DR/3y = 0 since )l has only a single component = H(x,t) 
by Eq. (3), and Vp = (5 by Eq. (1). 

Thus, Eqs. (1) - (3) lead to the following initial-boundary-value problem 

-V 

for the transverse velocity wave v(x,t) in the y-direction propagating in 
the x-direction, as a result of the sudden wall motion in the plane x = 0: 


3 v _311_ 

P o It dx 


(5) 

on n _ _ p 

01 t 'i 3 x 

♦ 

<&) 

v(x = U,t) = v H(t), t 1 0 

> 

(7) 

o 

A 

o 

1! 

/**N 

o 

II 

> 

> 

(8) 

> v (x, t = 0)/3t = 0, x • 0 


(9) 


where H(t) =0, L - - 0, and H(t) = i, t - + 0. 

Equations (5) - (6) represent a hyperbolic system, from which one obtains 
by elimination wave equations for the stress component 11 E n(x,t) and the 

*y 

velocity lield v(x,t): 


3 7 II 


= c ? - 
fix 


;v ri 

ot 7 


1 Oli 

t 3t 


( 10 ) 






whore 


C = (;,/p T) 2 (12) 

o 

is the (maximum) spued of the stress relaxation wave. Both i. (x,t) and v(x,t) 
satisfy similar (hyperbolic) wave equations with the same wave speed c. In 
the limit, t -> 0 and c -► ", with c ? t ► i /( , Eqs. (10) and (11) reduce to 

parabolic equations, according to which boundary values of Il(x,t) and v(x,t) 
would diffuse with infinite speed into the fluid (conventional Navier-Stok.es 


theory). Accordingly, only for t > 0 and 


■, transverse or shear waves 


exist in the fluid which represent, therefore, stress relaxation waves. 

According to Eq. (11) and Eqs. (7) - (9), the velocity field 
v(x,t) ~ v u(f,t) of the stress relaxation wave under consideration is des¬ 
cribed by the dimensionless initial-boundary-vnlue problem: 

3 2 u 3u _ 3 2 u 

+ 3* ” al 7 , (13) 

u(4 = 0, t) = H(t), t 0, , (14) 

u (f, t = 0) = 0, f, > 0 , (15) 

3u(C, £ = 0)/D t = 0, 4 > 0 , (16) 


where 

U(\,t) = v(x,t)/v, 4 = x/cr, t = t/t 

Equation (13) - (16) are solved by means of the Laplace transform 

i • 6) 

technique — winch gives 

u(4,s) -■ L(u(f.,f)] = / e s! " u(4,()d( 

0 

u((),s) = L[u(n,^) ] = J e H(Od( = s 


F/6 20/9 
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d 2 u 

w 


Since the initial conditions (15) - (16) vanish, Eqs. (13) - (14) yield for 
the transformed velocity u(£,s) the ordinary boundary-value problem, 

(s 2 + s)u = 0 , (20) 

u(£ = 0,s) = s ^ . (21) 

Since u(£,s) must be finite for £ -*■ “, the solution of Eqs. (20) - (21) is 

( 22 ) 

The inverse Laplace transform gives for the velocity field the complex 
integral ,. 

° 7 + 1 “ , 

_1 _ /e2 i 2 . e/* 

+ s; f ’e Sl ds. (23) 


- . . -1 -(s 2 + s)^ r 

u(£,s) = s e 


u(£,t) = ( 1 / 2 * 1 ) J s e ^ 

Y~i“ 

Hence, 

u (£,*) = 3$(£,-t)/3£ 

where 

y+ioo 

< Kf,,'t) = - 2 ~r / s' 1 f (£,s)e S *ds 


Y-i 00 


(24) 


(25) 


and 


f(£,s) = e (S2+S) W + a)* . (26) 

According to a known inversion integral,—^ the inverse transform of Eq. (26) is 

(27) 


fU,L) = r 1 lf(C,s)] = l o (h(t 2 - £ 2 )S H(t - £) 


where I U) is the modified Bessel function of order v. By Eq (25), 


<!>(£,£) = - L 1 [s 1 f(£,s)] = - J L X [f(£,s)]di 

t ° 

= - / f (£ ,-£)dC 
0 




i.e. 


♦(£,*) = - H u - o J e'^^I (^(a 2 - C 2 ) *)da 


(28) 


(29) 


From this potential, the dimensionless velocity field is obtained as 

u(£ ,t) = H(* - C) [e”^ * + hK J e~ Hu (a 2 - Z 2 )~ H I,(>s(o 2 - 5 2 )*‘)do] (30) 

£ 
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l 


in accordance with Eq. (24). By Eq. (17), the corresponding dimensional 
solution for the velocity field is 
v(x,u) = vH(cl - x) {e 

t/1 , , i 

+ (x/2ct) / e"' ?U [cc 2 - (x/ct) 2 ] 2 I^lu 2 - (x/ct) 2 ] 2 )daj. (31) 

x/cx 

Eq. (31) indicates that the transverse stress relaxation wave is dis¬ 
continuous at x = ct, *-he position of the wave front at time t. At any time 
0 i t S, only the region 0 1 x 1 ct of the fluid is excited by the wave, 
since v(x < ct, t) >0 and v(x > ct, t) = 0 by Eq. (31). The velocity signal 
v(0,t) = v H(t) generated at the boundary x = 0 at time t is thus transported 
with finite speed c = (p /p o t) 2 < °° in form of a discontinuous wave into the 
fluid space x > 0 as t > 0 increases. 


2 

Application of the asymptotic formula 1 (z) ~ e /(2:rz) 2 , |z| >> 1, and 
expansion of 2 = ^[a 2 - (x/ct) 2 ] 2 for large ;x-values in Eq. (3.1) yields, in 
the limit t -> 0, ct -* 0: 

00 2 1 

v(x,t) = v(2//tT) / e 6 dE, n = x/2(— t)' 2 . (32) 

n ‘o 

This is the familiar solution of the parabolic Navier-Stok.es equations due 

to Stokes. Eq. (32) suggests that v(x,t) > 0 throughout the entire fluid 

0 • x * for any, no matter how small time t > 0. Thus, the parabolic 

Stokes solution gives a completely misleading picture for a shear wave in 

form of a diffusion process which spreads with infinite speed. 

0 ] 2 

Fig. 1 shows u(c,t) versus £ for t = 10 , 10 , 10 , with wave fronts at 
t = 10 , 10 , 10~. It is seen how the perturbation u(0,i) = H(f) produced at the 
wall - 0 moves in form of a discontinuous wave into the fluid space Z - 0 
so that an increasing but finite region 0 1 (, - r, of the fluid is set into motion 

with increasing t. in tile limit ( = u(f.,0 = 1 throughout the fluid, 

< < 

j 


0 
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INITIAL-BOUNDARY-VALUE PROBLEM FOR v(x,0)-PULSE 


Another fundamental method for shear wave generation makes use of a velocity 


pulse v o =v(x,0)e^ generated of time t = 0 within a limited region |x| < Ax. 
The decay of this velocity pulse occurs in form of a shear wave with velocity 
field v = {0, v(x,t), 0} in the y-direction propagating in the x-directions. 

f -► -V •> •<—f 4 <■-> 

Accordingly, V• v = 3v/3y = 0, v-Vv = 0, and vVIl =v3Ii/Dy = 0, since II has 

only a single component II = n(x,t). Again, the stransverse wove "behaves" 

xy 

incompressible and linear, and Vp = 3 by Eq. (1). 

As in the previous problem, Eqs. (1) - (3) give the wave Eqs. (10) and 
(11) for n(x,t) and v(x,t), respectively. Hence, the shear wave produced by 
the velocity pulse v = v(x,0)e is described by the initial-boundary-value 

a y 


problem: 


3 2 v 1 3 v 2 
Ft 2 T Ft = c a? 


v ( x , t = 0) = V Q (x) 


I v(x, t = 0) /U t = w Q (*) , M 


< 

— CO 


(33) 

(35) 

(35) 


where |w (x)I - 0 is included for reasons of generality. The solution of 
o 

g) 

Eqs. (33) - (35) is accomplished by means of Riemann's method, — 

_ x+ct 
v (x, t) = e T (i;[v o (x-ct) + v^(x+ct)] + % / Y(x, t,a)da} (36) 

x-ct 

where 


T(x, t, ct) = v Q (a) (t/2t) lj_( 2 “ [c 2 t 2 - (u - x) 2 ] 2 ) / [c 2 t 2 - (a - x) 2 ] 2 
+ - [w (a) + tv- 1 - v (a)] I (^- [c 2 t 2 - (a - x) 2 ] 2 ) . (37) 

CO o o ZCT 

As a concrete example tor the initial conditions in Eqs. (34) - (35), an 
initial velocity distribution of the form of a Dirac pulse is chosen, 

v (x) = v 6(x), w (x) = 0, ™ . (38) 

o o o 
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In this case, the general solution in Eqs. (36) - (37) becomes, in dimen¬ 
sionless form, 

u U,t) = % e~ l [6(5- t)+ 6(5+ t) 

+ Ijilt?- e, 2 ] h +l Q ([t 2 - 5 2 ] 4 )], 151 1 

= 0 , Id ' *, (39) 

where 

u(5,t) = v(x,t)/(v o /2ci), 5 = x/2cx, t=t/2i . (40) 

Eq. (39) indicates that the shear wave spreads in the space |5| i ■» in 
form of a symmetrical wave, u(- 5,-£) = u(+ E,,t) due to the symmetry of the 
initial conditions (38). The wave is discontinuous at its fronts 5 = ± t, 
which propagate with the speed 

v(x,t) = ± c, x = ± ct . (41) 

In the limit c ■+ 0, ct -► 0, application of the asymptotic formula 

2 J, L 

I^(z) ~ e / (2nz) 2 , |z| >> 1, and expansion of z = [a 2 - 5 2 ] 2 for large '(-values 
in Eq. (39) yields 

u(5 »-0 = (2irf)- % 151 - - . (42) 

This is the corresponding solution of the parabolic Navier-Stokes equations. 

Eq. (42) would indicate that the shear wave has the form of a Gaussian ex¬ 
tending from 5 = _0 ° to 5= + °° for any, no matter how small time t > 0 (cor¬ 
responding to an infinite speed of propagation). It is obvious that the 
solution (42) is physically not meaningful. 

In Fig. 2, the dimensionless velocity field u(5,-t) of the shear wave is 
shown versus 5 for t = 10®, 10^, 10 2 , the wave fronts being in each case at 
5=1 t. Due to the finite wave speed c, the fluid is not excited in the 
region |f| > f ahead of the wave fronts. The shape of the wave is flat 


% 
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with relatively steep flanks leading to the discontinuous fronts. Thus, 
the shear wave does not resemble the Gaussian of the parabolic theory, Eq. 


i 

1 


i 


(40). 
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CONCLUSIONS 

A generalization of the Navier-Stokes equations is presented considering 
viscous stress relaxation, which results in a physically meaningful theory for 
transverse waves in viscous fluids. The fundamental speed of the stress 
relaxation waves is given by c = (p/px) J , where \i is the viscosity, p is the 

9) 

density, and t is the relaxation time of the stress tensor. For any medium, 

it is c ~c , where c is the speed of the longitudinal waves, e.g., 

c = 1.3 x 10cm sec ^ and c * 1.5 x 10^ cm sec ^ for water at T = 20°C and 

s 

P = 1 atm. 
o 

Exact solutions are derived for stress relaxation waves propagating in 
the x-direction due to velocity pulses v(0,t) and v(x,0) in the y-directlon, 
respectively. For the geometry of these transverse waves, the nonlinear, 
generalized Navier-Stokes equations become linear, so that the solutions given 
hold for waves of arbitrary intensity. The solutions are discontinuous at 
the wave fronts, which is typical for hyperbolic field equations. The corres¬ 
ponding solutions of the conventional Navier-Stokes equations indicate a 
diffusion process with infinite wave speed, i.e., give a qualitatively and 
quantitatively insufficient picture of the propagation of transverse waves 
in fluids. 

In the simplified stress relaxation equation (3) proposed, the term 

'[•Vv neglected since it is of the order-of-magnitude of (h/t)jVv| 2 , 

which is nonlinear in the derivatives. It should be noted that the term 

H*Vv vanishes exactly for the wave problems treated above, Il*Vv=0 since v = 

{0,v(x,t),0} and II has only a single component n (x,t). For this 

xy 

reason, the solutions presented arc exact solutions of the nonlinear 
Navier-Stokes equations with viscous stress relaxation. 
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